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Statistical Distribution Patterns of 
Ocean Waves and of Wave-Induced Ship 
Stresses and Motions, with 
Engineering Applications 


By N. H. Jasper,? MEMBER 


Little is known about the frequency of 
occurrence of the various magnitudes of 
ocean waves. Even less is known about 
the severity of the wave-induced motions 
and stresses which ships experience in 
service. The intent of this research effort 
is to show that, by utilization of statistical 
methods, it is possible to describe and 
predict service conditions for ships in an 
orderly and relatively simple manner de- 
spite the general complexities of the prob- 
lem. Wave observations taken continu- 
ously over a period of 6 years at several 
weather stations in the Atlantic Ocean 
were studied. Wave-induced motions and 
stresses in ships obtained under a wide 
range of operating conditions were studied 


INTRODUCTION 


It is probably correct to state that the ship- 
building industry knows less about the service 
conditions under which its product must operate 
than does any other major construction industry. 
In contrast to the aircraft builders, shipbuilders 
have made only a small effort to establish actual 
service stresses and motions or in incorporating 


1 The paper given here was originally submitted to the School of 
Engineering and Architecture of the Catholic University of America, 
Washington, D. C., in partial fulfillment of the requirements for the 
Degree of Doctor of Engineering. 

2 Engineer and Deputy Division Head, Vibrations Division, 
David Taylor Model Basin, Navy Department, Washington, D. C. 

Presented at the Annual Meeting, New York, N. Y., November 
15-16, 1956, of THe Socrtety oF NAVAL ARCHITECTS AND MARINE 
ENGINEERS. 


for seven different ships. On the basis of 
an analysis of voluminous experimental data 
it is concluded that the probability dis- 
tributions of wave height, wave length, 
wave-induced pitch, roll, and heave mo- 
tions of ships, and wave-induced ship 
stresses may all be approximated by a one- 
parameter-type distribution function when 
the environmental conditions are steady, 
whereas these variables will tend to follow 
the two-parameter logarithmically normal 
distribution when the environmental con- 
ditions are allowed to vary over a wide 
range. Applications of the results of this 
study to the solution of engineering prob- 
lems are described with the help of perti- 
nent illustrations. 


the results of such studies into methods of ship 
design. A deterrent to progress in this direction 
has been the complexity of the problem as evi- 
denced by the extensive variability of the sea 
and of the response of the ship to the sea. The 
intent of this paper is to show that, by utilizing 
statistical methods, it is possible to describe and 
predict service conditions for ships in an orderly 
and relatively simple manner, even though the 
problem appears formidable. 

The experimental approach taken in this paper 
should complement the theoretical work of St. 
Denis and Pierson although the present study 
was developed independently of theirs and is not 


TABLE 1 STATISTICAL TERMS AND DEFINITIONS 


It should be noted that in the present work only estimates of the statistical parameters can be 


given 

variate Property under study, for example, magnitude of stress variation, height of 
waves. 

N 
Ze, 

E The mean square value of z, E =4=1— 

P(%a) Probability of the event « < xa with respect to the random experiment «. The 
probability P means that in a long series of repetitions of ¢ it is practically 
certain that the relative frequency of occurrence of x = Xa is equal to P where 
0 <P <\J. A value of P = 1 indicates that the event will occur with every 
experiment e, P = 0 indicates that the event will never occur. 

p(2) =4 p(z) 

dz 

population Aggregate of all the members of the distribution subject to a statistical 
study. The members may be individual measurements or the results of indi- 
vidual experiments. 

sample Number of N members drawn (usually at random) from the population. Here 
N is the ‘‘sample size.’’ 

c.d.f. Abbreviation of the term ‘‘cumulative distribution function’’ P(x) which 
denotes the probability that the variate z will take on a value less than or 
equal to 2, i.e., P(z,) =Ple<-2,). 

d.f. Abbreviation of the term ‘‘distribution function’’ p(z) which is the derivative 
of P(x) with respect to g, i.e., 

+ 00 
p(2) = P(z) and { p(z) de = 1 
dz 

v = oe 6 

EA 
N 
2 2, 

mean value of z ee where z, is the individual member of the population or sample and N 
is the number of members of the population or sample. 

z Mean value of a sample of values of z 

é Mean value of the population. 

: i (x; - 2)? 

variance of z = aia a a The variance is a measure of the dispersion of the values of 
z. The greater the variance the greater the dispersion. 

a? Variance of the population. 

standard deviation Square root of the variance 


quantile or fractile Corresponding to a given value of the c.d.f. P(x), there is a value of the 
variate below which a fraction P(z) of all the members of the population lie. 
The value of this variate is the quantile or fractile corresponding to P(z), 


confidence bands Interval within which the ‘‘true’’ distribution will fall with a certain probability. 


random A method of drawing a sample when each member of the population has an equal 
chance of selection. 


limited by the assumptions required by their 
mathematical analysis. 

The purpose of this paper is to present evidence 
that the wave-induced pitching, rolling, and 
heaving motion of ships, as well as the associated 
hull-girder stresses follow simple, mathematically 
defined statistical frequency or probability dis- 
tributions. It will be shown that the distribution 
functions of motions and stresses are of the same 
type as those applicable to the height of ocean 
waves.’ In particular it will be shown that the 
distribution function approximates the single 
parameter “Rayleigh” type for a given state of 
the sea, ship speed, and course. On the other 
hand, when distributions of wave height, ship 
motion or stress experienced over a wide range 
of conditions are considered, the two-parameter 
“Jogarithmically normal’’ distribution is shown 
to be applicable over the range of conditions 
observed thus far. The latter distribution is the 
result of a summation of many distributions of 
the first type. Finally, the extreme values of 
ship motions and stresses will be examined to 
evaluate the feasibility of the prediction of ex- 
treme values. 

For the sake of brevity, the distribution ap- 
plicable to uniform conditions of the sea, ship 
speed, and course will be called ‘‘short-term”’ 
distribution, whereas the function which repre- 
sents the distribution obtained when the seas, 
ship speeds, and courses are allowed to vary over 
a range of conditions, will be designated as a 
“long-term’’ distribution. 

The conclusions arrived at in this paper are 
based on numerous observations of wave height, 
hull-girder stress‘ (for seven ships) and ship 
motions (for five ships). The experimental data 
were obtained over a period of years and cover 
a wide range of service conditions. 

There are many applications in which a knowl- 
edge of the frequency distribution of hull mo- 
tions, stresses, and the heights of ocean waves can 
be used to advantage. For example: 


1 Prediction of the most probable amplitudes 
of roll and pitch motion of ships under given 
environmental conditions. 

2 Estimation of the extreme values of ship 
response or of wave heights encountered over 
given periods of time. 

3 Statistical estimation of the capacity for 
which shipboard stabilization equipment must be 
designed. 


3 The study of ocean waves originated as a by-product of the basic 
investigation of the ship’s response to the sea. 

4 Wherever the term stress is used in this paper, a stress computed 
from the measured strain is to be inferred. 


4 Estimation of the endurance strength of the 
ship structure. 

The experimental work carried out in connec- 
tion with this problem has been a co-operative 
effort of many people and several organizations. 
In particular the U. S. Navy, The Society of 
Naval Architects and Marine Engineers, the, 
U. S. Coast Guard, the U. S. Weather Bureau, 
and the Esso Company have contributed much 
to make this research possible. The special 
instrumentation and methods used to collect the 
data, as well as some of the sea tests are described 
in references (1-4). This paper will not be con- 
cerned with the methods with which the data were 
obtained except in so far as they are pertinent 
to the statistical analysis. 

The general plan of presentation will be as 
follows: Some pertinent statistical background 
information will be given first, followed by a short 
outline of the sources of the data employed in this 
study. The analyses of short-term and long- 
term distributions will be given separately. Each 
analysis applies statistical tests of significance to 
the hypothesis that the experimental data are 
samples from a specified distribution. The re- 
sults of the analyses are discussed and applica- 
tions of the results to practical problems are 
illustrated, followed by a summary of the major 
conclusions. 


STATISTICAL BACKGROUND 


The wave heights, the variations in ship mo- 
tions, and hull-girder stresses experienced under 
a given set of conditions will be described in terms 
of their distribution functions. The statistical 
methods applied in the present study comprise 
techniques which were used to obtain, classify, 
analyze, and present large masses of experimental 
data. In particular they were used to: 

1 Collect the data by sampling procedures. 

2 Classify and condense the data. The data 
were classified according to the magnitude and 
frequency of the measured variations. Each vari- 
ation is understood to mean a fluctuation in the 
magnitude of the variable from the largest posi- 
tive value to the succeeding largest negative 
value; for example in the case of a rolling motion 
the variation would denote a motion of the ship 
from the extreme starboard position to the ex- 
treme port position. 

3 Present the data. The data are given in 
tabular form as well as in the form of statistical 
distribution patterns, such as histograms and dis- 
tribution functions. 

4 Analyze the data. With the data presented 


5’ Numbers in parentheses refer to the Bibliography at the end 
of the paper. 


Gage Locations 


ee 


LBP 528’ Loaded Displacement, Gages: 6” below main 
Beam 76° 21,050 Tons deck on side plating; 
Depth 44%’ 89’ fwd. Frame 106. 


SS GOPHER MARINER, C-4, Dry Cargo 


Gage Locations 


al Se ee 


LBP 503’ Loaded Displacement, Gages; on main deck, 
Beam 687 21,880 Tons 9%,’ inside rail, at 
Depth 3914” Frame 57. 


SS ESSO ASHEVILLE, T-2, Tanker 


Gage Location 09 


be Lg 


Lis 67 Loaded Displacement, Gage; approx., upper 
Beam 56°10%* 13,000 Tons deck, 12’ stbd. of 
Depth 37°74” centerline, Frame 92. 


SS OCEAN VULCAN 


eP=F) 


Gage Location +--+, 


Gage Location Main Deck 


LBP 306” Loaded Displacement Gage; on centerline, 

Beam 367 1630 Tons main deck, 18°’ fwd. 

Depth 20.67° Frame 90. 10.717 
USS FESSENDEN, DER-142 above neutral axis. 


Longit. Hull Girder Stress 
at Amidships 


Heave Acceleration 


Loaded Displacement, 
As OE 2,500 Tons 
LBP 300° 
Beam 4173/4’ 
Depth 24° 74”" AP 


USCGC UNIMAK 


Fic. 1 LocaTIon oF STRAIN GAGES 
(Circles indicate gage locations.) 


in a suitable manner, analytic functions are fitted 
to the data, and tests of significance are applied 
to determine whether certain statistical predic- 
tions and conclusions can be reached with con- 
fidence on the basis of the available data. 

Some basic terms used in statistical work are 
defined in Table 1. A distribution function indi- 
cates the relative frequency or probability of oc- 
currence of a particular event. The best known 
type of distribution is the “normal” distribution. 
It is theoretically possible, but not often prac- 
ticable, to convert distributions which are not 
“normal” into a normal distribution by a change 
of variables. One such conversion which will be 
utilized here is to take the logarithm of the vari- 
able. Thus the “log-normal”’ distribution is ob- 
tained, which signifies that the logarithm of the 
variable is normally distributed. Commercial 
chart paper (probability chart) is available which 
is designed so that the integral of the normal or 
log-normal distribution will plot as a straight line 
on the chart. The log-normal law has been 
found applicable in various diversified fields. 
See, for example, the extensive bibliography given 
by V. T. Chow (5). 

The probability P(x) takes on definite values for 
a given population. An estimate of P(x) is often 
desired on the basis of samples of the members 
of the population. Statistics provide the tools for 
making such estimates. It will be assumed that 
the sample is taken from a population for which 
the quantity under discussion is defined by a 
certain analytical expression. The reliability 
of this hypothesis is evaluated by the application 
of statistical tests of significance which enable one 
to estimate the expected variation of the measured 
data from the analytically defined values, if the 
hypothesis is true. 

In order to determine whether the variation of 
the sample data from the assumed distribution is 
consistent with the hypothesis, use will be made 
of a theorem in statistics (6) which states that 
the sample quantile® x, is asymptotically normal 
with mean value £ and standard deviation 


I ea = 2) 
p N 


where NV is the size of the sample, p the probability 
density, P the cumulative probability at x = xp 
and & is the corresponding quantile of the popu- 
lation from which the sample is drawn. Other- 
wise stated: If a sample of N-values is chosen 
at random and if these values are numbered in the 
order of increasing magnitude, then the magni- 


_ ° A quantile denotes the variate corresponding to a given fraction; 
i.e., the 0.40 quantile is that value of the variate below which 40 
per cent of all values lie. 


tude of any particular value, say the rth cannot 
be predicted in advance of the drawing of the 
sample; there exists a certain probability function. 
If N is large, and if neither r/N nor (N — r)/N is 
too small, then the probability function for the rth 
value is approximately normal in form. If JN is 
increased and if 7 is also increased so as to keep 
r/N as nearly constant as possible (as for a fixed 
quantile of the sample) then the approximation 
to normality increases, becoming perfect in the 
limit as NV —~ o; the width of the probability 
spread simultaneously decreases, becoming zero 
in the limit. 

If a curve is plotted on each side of the analytic 
cumulative distribution curve at a distance of 
1.65 standard deviation of the distribution of the 
corresponding quantiles, then on the average nine 
tenths of the observed quantiles would be ex- 
pected to fall between these limits and 10 per cent 
of the measured quantiles would be expected to 
fall outside these limits. This control or confi- 
dence curve is therefore designated as the 90 
per cent significance level’ and indicates the 
probability of a random value falling within the 
limits when the hypothesis is true. In order for 
the deviation between the hypothetical distribu- 
tion and the sample values to be considered not 
significant at this level, it would be necessary to 
have less than 10 per cent of a large number of 
observed values fall outside these limits. In 
practice it has been found that wrong conclusions 
will seldom be drawn if the level of significance is 
set at 90 per cent or less; that is, the hypothesis is 
probably true if the deviations are no greater than 
could be expected, on the average, in one out of 
10 tests. On the other hand if the deviations lie 
outside the 99 per cent limits, that is, if they are 
larger than would be expected once in a hundred 
similar tests, then the hypothesis should be re- 
jected confidently. A distribution may be of 
practical significance even though the statistical 
tests indicate a significant deviation between it 
and the measured data. 

In many cases a visual inspection of the scatter 
is sufficient to indicate that the assumed distri- 
bution is an acceptable representation of the actual 
distribution. The types of basic distributions 
that will be considered here are the normal, log- 
normal, and the Rayleigh distributions. 

The normal distribution of x is defined as fol- 
lows: 

1 =(e=5? 

= 3 AP — @ SUS ap @ 
a V 20 


p(x) = 


7 Some statisticians prefer to label the limits according to the 
probability with which a random measurement would fall without 
the limits; i.e., the 90 per cent significance limit would then be re- 
ferred to as the 10 per cent level. 


P(x) = Spx) dx 


where 


p(x) is the probability density of x 

x is the variate 

E is the mean value of the variate 

o? is the variance of the population 
Thus the two parameters £ and o define the dis- 
tribution completely. 

The logarithmically normal or log-normal dts- 
tribution is a normal distribution of the logarithm 
of x, and 

TL = (log 2 = u)? 
De) xo WV 20 : 
where u is the mean value of log x and a is the 
standard deviation of log x. 
The Rayleigh distribution is defined as follows 


p(x) = Be-*/F, x 5 0 


where E is the mean value of the squares of x; thus 


TasLE 2 USCGC Unimak—CompuTaTION oF CONFI- 
DENCE LIMITS FOR DATA SHOWN IN FIGs. 2 AND 3. 


(Pitch angle, sea state 5, head seas, ship speed 71/2 knots.) 


z P P oa 

Variation in | Probability Density | Probability | Standard Beviation eG Zee 
Pitch Angle (per degree) z of Fractile Upper Limit | Lower Limit 

deg ‘a 2 Pal pde “Der | ee 
rg oie ltd ns) N 67% Confidence Limits, deg 

0.5 0.0210 0.050 0.2186 0.72 0.28 

1 0.0413 0.0206 0.2239 1.22 0.78 

2 0.0775 0.0809 0.2290 2.23 1.77 

3 0.1046 0.1729 0.2353 3.24 2.76 

4 0.1193 0.2923 0.2482 4.25 3.75 

5 0.1244 0.4097 0.2573 5.26 4.74 

6 0.1190 0.5294 0.2730 6.27 5.73 

7 0.1051 0.6441 0.2965 7.30 6.70 

8 0.0875 0.7407 0.3260 8.33 7.67 

9 0.0688 0.8188 0.3644 9.36 8.63 

10 0.0511 0.8786 0.4160 10.42 9,58 

ll 0.0362 0.9221 0.4820 11.48 10.52 

12 0.0243 0.9519 0.5732 12.57 11.53 

14 0.0118 0.9784 0.8019 14.80 13.20 

16 L 0.0052 0.9911 | 1.1756 17.18 14.82 

The 67 percent confidence limits are given by x+ g; the 90 percent by x + 1.6Sa, the 95 percent 
by x + 1.96; and the 99 percent by x + 2.58.0. 

The value of the standard deviation computed in the table assumes that there are no errors intro 
duced in the process of measurement and reading of recorded data. The variance of the latter errors 
in roughly @2 = 0.30. Thus the confidence limits would be broader than given In the table if an 
allowance for these errors is made. The following formula may be used: 

corrected variance = 0 2= 07 +7 | 


Test Conditions: Sea State 5, Head Seas, Ship Speed 7 !2 knots 
Significant Wave Height 21 ft. 
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percent/degree pitch angle 


p= 


ee ab 


The ship pitched at a rate of about 330 variations per hour. | 
The test sample consisted of 236 variations in pitch angle. 
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Fic. 2 USCGC Unmak, DIstRIBUTION OF VARIATION IN PITCH 
ANGLE—SAMPLE 1 
(This sample is representative of 129 samples that were analyzed.) 


the single parameter E defines the distribution. 
It is possible to represent the distributions cor- 
responding to different values of E by a single 
distribution, employing a change of variables. 
Thus let 

Gs 


VE then ./E p(x) = 2uve7” 


v= 


Integrating to obtain the cumulative probability 
VE Si" pe) dx = VE P(e) = VE lL —e-*) 
or 

IAGy) Sl =e] 


and 
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v2 = 2/E 


Fic. 3 USCGC Unimaxk, CuMULATIVE DISTRIBUTION OF VARIATION IN PitcH ANGLE—SAMPLE 1 
(See Fig. 2 for corresponding histogram. Note: The plotted confidence limits do not include an allowance for errors of measurement.) 


l= JAG) = Ge 


Therefore if the variable under study is expressed 
in terms of the standardized variate v then the 
function [1 — e—°?] will represent the cumulative 
distribution of all Rayleigh distributions. This 
cumulative-distribution function will plot as a 
straight line if the logarithm of [1 — P(x,)], 
that is, the probability of exceeding x,, is plotted 
against v? as has been done in Fig. 3. 

It should be noted that the methods of sta- 
tistics are strictly valid only if the measured val- 
ues are independent samples. The physical 


5 The USCGC Unimak, a former U.S. Navy 
AVP class vessel, now fitted as a weather ship, 
operating at weather station B in the North 
Atlantic Ocean. Some data also were obtained 
for a sister ship, the USCGC Casco. These 
vessels are approximately 300 ft long and displace 
about 2500 tons. 

6 A destroyer approximately 380 ft long oper- 
ating in the Atlantic Ocean. 

7 An aircraft carrier approximately 900 ft 
long operating in the Atlantic Ocean. 

The actual measurements of strain obtained 


p = Percent/degree Pitch Ancle 


The ship pitched at a rate of about 370 variations per hour. 
The test sample consisted of 175 variations in pitch angle. 


T T AV 1 T T 
Test Conditions: Sea State 5, Quarter Following Seas, 
Ship Speed 7'2 knots, Significant | 
Wave Height 21 ft. 


Theoretical Rayleigh Distribution | 
E = 20.2 (degree)¢ 


Experimental Histogram 
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Fic. 4 USCGC Unimax, DistRIBUTION OF VARIATION IN PitcH ANGLE—SAMPLE 2 


(This sample represents the poorest fit to the Rees distribution of 129 samples that were an- 
alyzed. 


data are not, in general, independent of one an- 
other. However if the tests of significance indi- 
cate, on the basis of assumed independence of all 
observations, that these data reasonably can be 
assumed to derive from the distribution speci- 
fied by the hypothesis, then this conclusion is 
even more sound than the tests of significance 
indicate. 
SOURCES OF DaTA 

Stresses 

Service stresses were obtained for the following 
ships: 

1 The SS Esso Asheville, a T-2 Maritime Com- 
mission design tanker, which operated in coast- 
wise service between the United States, Central 
and South America. 

2 The SS Gopher Mariner, a C-4 Maritime 
Commission design dry-cargo vessel, which oper- 
ated in the transatlantic service. 

3 The SS Ocean ‘Vulcan (7), a dry-cargo vessel. 
No specific data are available to the author as to 
its area of operation. 

4 The USS Fessenden, a destroyer escort, oper- 
ating in coastal waters off the U.S. East Coast. 


for these ships, excepting the Ocean Vulcan (7), 
were made by the U.S. Navy’s David Taylor 
Model Basin. In many cases these measure- 
ments were obtained by means of automatic 
instrumentation which did not require the pres- 
ence of an operator aboard the ship. The strain 
gages were installed on the strength deck, near 
amidships, see Fig. 1. The tests of six of these 
ships (Gopher Mariner, Esso Asheville, Unimak, 
Fessenden, Destroyer, Carrier) covered a period 
of about 3 to 8 months, predominantly in the 
winter season; the Ocean Vulcan data however 
cover a period of 214 years. 


Motions 


Ship motions were obtained by Taylor Model 
Basin personnel on the Esso Asheville, Casco, 
Unimak, a destroyer, and an aircraft carrier. 
The motions were recorded by means of an auto- 
matic recorder which took a record at preset 
intervals of time for a predetermined duration of 
time. The complete tests on the Esso Asheville, 
Casco, and Unimak are described in references 
(3, 1 and 4), respectively. 
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(See Fig. 4 for corresponding histogram.) 


1-P Percent 


Wave Data 


Observation of significant wave height® and 
wave period, obtained over a period of 6 years, 
have been made available by the U.S. Weather 
Bureau. These observations were made every 
3 hr by trained weather observers in accordance 
with instructions prescribed by the U.S. Weather 
Bureau (8). Only one quantitative measure- 
ment was reported each time the sea was ob- 
served. Darbyshire (9) obtained measurements 
of wave heights, at Ocean Stations I and J by 
means of a wave meter which has been developed 
in England. Thus one set of measurements is 
available which may be used as a check of wave- 
height observations reported by weather ships. 


ANALYSIS OF SHORT-TERM DISTRIBUTION 
FUNCTIONS 


It is the purpose of this section to show that, 
on the basis of the experimental evidehice, wave- 
induced motions and stresses of ships are dis- 
tributed according to the Rayleigh distribution 
with a high degree of statistical significance pro- 
vided the sea conditions, ships speed,.and head- 
ing remain uniform throughout the period under 
consideration. 

It has been accepted generally that the heights 
of ocean waves passing any given point follow the 
Rayleigh distribution, provided the associated 
power spectrum has most of the energy con- 
centrated in a narrow frequency band. This dis- 
tribution is valid if measurements of wave height 
are taken over a relatively short period of time, 
of the order of 1 hr, during which interval the sea 
conditions do not change appreciably. It can 
be shown that this distribution is the same as 
that representing all wave heights, in the area 
under consideration, at one instant of time. The 
statistics of the Rayleigh distribution were applied 
by Longuet-Higgins (10) to the analysis of the 
distribution of wave heights in a sea character- 
ized by a narrow spectrum. This work was ex- 
tensively. applied by Pierson, Neumann, and 
James (11) to the study of ocean waves. 

Inasmuch as hull stresses and motions are 
caused by the passage of the waves past the ship 
it will be appropriate to test the hypothesis that 
the responses of the ship to the waves, i.e., stresses 
and motions, may similarly be statistically repre- 
sented by a Rayleigh distribution. 

For this purpose a number of tests? were run 


8 Significant wave height is the average height of the upper third 
highest waves. 


9 Each test yielded a continuous !/2-hr oscillogram record of ship 
motions and stresses applicable to a specific combination of condi- 
tions. 
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TABLE 3 CHI-SQUARE TEST APPLIED TO TEST GOODNESS 
OF Fit OF RAYLEIGH DISTRIBUTION TO TEST RESULTS 


Expected 


Number of 
Class Limits | Variations | Number of 
Measured | Variations 


UNIMAK 0- 3 deg 40 48 
eee | eo | a la 
6 8 49 50 
8-10 33 33 
10-14 a1 26 0.04 
14-19 5 2 45 
0- 21/8deg| 46 35 3.46 | Ex? = 117 
21/8- 31/8 25 33 1:93) | alrihit 
31/8: 41/8 31 34 0.26 
41/8 51/8 rE) 28 0.89 
51/8 61/8 14 20 1.80 
61/8- 81/8 21 19 3.36 
81/8-101/8 6 6 0 


Value of 


Remarks 
Chi Square 


27 =6.4 
Good Fit 
P=0.15 


Sample Tested 


UNIMAK 
Pitch Angle 


See Figures 4, 5 P= 0,92 


Ey? = 4.4 
Very Good Fit 
P= 0.30 


DESTROYER 


Pitch Angle 1-2 34 41 1.20 
See Figures 6, 7 ; 


7-10 11 15 1.07 
AIRCRAFT CARRIER 0- 5 deg 43 45.96 0.19 
Pitch Angle 9 0.59 
See Figures 8, 9 Sao ue ee) # 

9-12 36 43.69 1.35 
12-15 32 27.93 0.59 
15-18 15 14.23 0.04 
18-23 8 7.42 0.05 
28.75 1.36 | Ex? = 6.8 


Good Fit 
0.00 | P=0.75 


Ly7=2.8 
Very Good Fit 
P= 0.30 


AIRCRAFT CARRIER 
Stress in Main Deck 
See Figures 10, 11 


0- 750 psi 35 

750-1500 68 66.47 

1500-2250 61 66.40 0.44 

2250-3000 31 41.29 2.56 

3000-3750 25 18.97 1.92 

3750-5500 10 8.13 0.8 
> 


on the USCGC Unimak under a wide variety of 
combinations of sea, speed and ships heading to 
the waves. The hull stresses, pitching, rolling, 
and heaving motions were measured for 39 sets of 
these combinations. The parameter E which 
defines the Rayleigh distribution 


p(x) = Oe g—/i 


was computed for each of the frequency distri- 
butions of motions and stresses measured under 
the various combinations of environmental con- 
ditions. The individual histograms were plotted 
and compared with the Rayleigh distributions 
corresponding to the computed value of E. A 
typical distribution, thus determined, is shown 
in Figs. 2 and 3. The pattern which evidenced 
the poorest agreement between the experi- 
mental and the Rayleigh distribution is shown 
in Figs. 4and 5. 

It has been shown in the section ‘Statistical 
Background,” that all Rayleigh distributions may 
be represented by a single analytic expression if a 
new variable v is used rather than x directly. The 
cumulative distribution function then plots as 


TaBLE 4 USCGC Unimax—PircH ANGLE—COMPARISON OF PREDICTED 
AND MEASURED MaxIMUM VALUES 


lr Sign Wave Heading of Ship | Sampling | Total No of | Computed | Measured Max| Predicted Max Procedure for Prediction of 
ae Height | Waves Relative] Speed} Time | Variations in| Value of E, Value Value Probable Maximum Value 
tate ft to Ship hate aia Sample, N (deg) | dep deg from a Rayleigh Distribution 
2 6 Head Seas 10 30 240 4.00 4.8 4.7 
14 32 296 4.48 5.3 5.1 Let X max = aust ane maximum 
Quarter 0 | 37 249 1.97 2.9 33 Saas SAR 7 
GEESE |! 1g 32 274 1.86 3.3 3.2 values of x 
Beam Seas | 10 29% 276 5.40 6.3 5.5 Then according to Longuet-Higgins 
14 32 343 3.75 5.1 4.7 
i x, = E% gt 
Quarts 14 29 282 1.68 3.9 31 max 
0 own where @= log. N — log. {i - a Qa = =| 
Polo: |! 74 32 308 1.86 3.8 3.3 
IL Seas IL 
3 7-9 Head Seas 10 26 253 10.35 7.3 7.6 Sample Size Xmax 
14 32 321 9.62 7.3 7.5 N e% 
17 27 277 9.59 6.3 7.3 |e’ [om |OS07Ia| 
Quarter 10 | 27% 269 1.79 31 3.2 a» | sce 
Head Seas 5 
17 37 36! 0 H H 
) 2.02 3.9 3.5 5 1.366 
ass Seas | o | 2 262 3.52 4.3 44 10 | 1.583 
varter 
Following | 17 31 270 1.85 3.9 3.2 2 sede 
Seas 50 2.010 
OS | iy |] 193 1.28 24 2.6 Ma | Sete 
| Seas | 200 2.323 
| 16 Head Seas IL rT] 30 289 26.12 10.3 12.2 500 | 2.509 
5 21 Head Seas 1% 21h 236 47.32 18.0 16.0 1000 2.642 
14 24 235 56.00 15.5 17.5 2000 2.769 
IL ee 
Quarter Th) 38% 347 35.13 13.5 14.3 5000 2.929 
Head/Seas) ||| 14 28% 274 47.14 16.5 16.6 10,000 | 3.044 
Beam Seas 7% 30 182 16.78 8.8 9.3 | 20,000 3.155 
Quarter 7% 2842 175 20.24 9.8 10.4 50,000 3.296 
earn ne 10 29 143 28.00 15 11.8 100,000 | 3.400 
14 27 136 17.2 8.3 9.2 
Following 1% 33 160 31.12 12.5 12.6 
S225 4 | 3i 140 20.54 10.8 10.1 
25¢ T =F T 7 \ =o Ley 


p = Percent/degree Pitch Angle 


—— 


The ship pitched at a rate of about 400 variations per hour. 
The test sample consisted of 270 variations in pitch angle. 


———— 


Test Conditions 
Significant Wave Height, 14 ft 
Head Seas 
Ship Speed, 15 knots 


Theoretical Rayleigh Distribution 
E=17.1 (degree)2 


Experimental Histogram 


0 1 2 3 4 5 


6 7 8 9 10 11 


Variation in Pitch Angle in degrees 
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a straight line on the type of probability chart 
that was devised in the section mentioned (see 
Fig. 3). 

A statistical test of significance was applied 
to the hypothesis: ‘The variations of the meas- 
ured sample from the assumed population, 
specified by E, are no larger than may be ex- 
pected according to the laws of chance.” 

Aside from the deviations attributable to chance 
it was necessary to consider the expected vari- 


11 


DISTRIBUTION OF VARIATION IN PitcH ANGLE, DESTROYER 


ations due to the inaccuracies of the measurement 
process. The variances of the fractiles for the 
sample under consideration were computed ac- 
cording to the theorem stated in the section cited. 
Fig. 3 which represents the results for a typical 
test indicates that the measured data fall well 
within the 90 per cent limits; a second set of curves, 
the 67 per cent limits also was computed and the 
data satisfy this more restrictive test. Thus a 
high degree of confidence can be placed in the hy- 
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Test Conditions: Significant Wave Height, 16 ft 
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40 The ship pitched at a rate of about 320 variations per hour. 
The test sample consisted of 210 variations in pitch angle. 
Theoretical Rayleigh Distribution 
E=6.3 (degree)? 
30 
Experimental Histogram 
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TABLE 5 STATISTICS OF LOG-NORMAL DISTRIBUTION OF WAVE- 
HEIGHT OBSERVATIONS IN THE NORTH ATLANTIC OCEAN 


h © Mean Value of Variance (02) of 
Period of | Ocean} | stitude Longitude Number of | Median Value of Logarithm of Logarithm of 
Records | Station aaa crn Observations® | Sign. Wave Height Sign. Wave Height| Sign. Wave Height 

1/49- a] A | 62°00°N | 33°00°W 12,891 6.34 1.8472 0.4524 
1/49-12/54 B 56°39°N | 5100 W 15,547 6.59 1.8859 0.4434 
1/49-12/54 c 52°45°N | 35°30W 16,857 6.75 1.9100 0.3763 
1/49-12/54 it) 44°00°N | 41°00°W 16,804 6.26 1.8338 0.3843 
1/49-12/54 E 35°00°N | 48°00W 16,777 4.56 1.5164 0.3765 
1/47- 6/53 | 61°00°N | 15°20°W 11,274 7.36 1.9960 0.3747 
1/47- 6/53 J 52°30°N | 20°00°W 12,016 7.40 2.0021 0.3863 
1/49-12/53 K 45°00°N | 16°00°W 11,182 6.20 1.8240 0.3033 
1/49-12/53 M 66°00°N | 02°00°E 14,324 4.99 1.6076 0.2344 
1/49- 6/54 H 36°00°N | 70°00°W 14,607 5.08 1.6251 0.4237 
*Each observation characterizes a sea state. 
All logarithms are natural logarithms. 


pothesis that this sample was taken from a Ray- 
leigh distribution. Consideration of the vari- 
ances due to instrument errors would have re- 
sulted in wider confidence limits, Table 2. 

As a further check on the validity of the hy- 
pothesis, the chi-square test was applied to the 
grouped data, in accordance with the method of 
chapter 23, reference (12). This test gave a 
rather low value of chi-square, Table 3, thus indi- 
cating that there are no inconsistencies between 
the test data and the hypothesis. 

Figs. 4 and 5 show the distribution for the poor- 
est set of data obtained in the sea tests of the 
Unimak. In this case also there are no signifi- 


cant differences between test data and the 
hypothetical Rayleigh distribution on the basis 
of either the chi-square test or the variation in the 
distribution of quantiles. 

Inasmuch as none of the more than one hundred 
distributions of ship motion or hull stresses (longi- 
tudinal, hull-girder bending moments) indicated 
significant deviations from the hypothesis, this 
hypothesis may be accepted confidently for ships 
of the type tested. 

Thus far only the results of the Unimak sea 
trials have been discussed. A similar analysis 
also has been made of the results obtained from 
sea tests (in the Atlantic Ocean) on a destroyer 
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(See Fig. 8 for corresponding histogram.) 
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and a large aircraft carrier. Figs. 6 through 11 
show typical distributions obtained from the sea 
trials of the destroyer and carrier. It is again 
evident that there is no significant difference be- 
tween the measured data and the hypothesis that 
the data are samples taken from a Rayleigh dis- 
tribution. 

It will be of interest to compare the largest 
measured value out of a given number of meas- 
urements with the most probable largest value 
predicted by the statistical theory. Longuet- 
Higgins shows that the largest probable value 
Xmax OUt of N-measurements is /E times a con- 
stant if the population is of the Rayleigh type 
where the constant is a function of NV only. For 
N large this constant becomes approximately 


equal to V log. N. Table 4 gives the comparison 
of predicted and measured values for a number of 
tests. There appears to be satisfactory agree- 
ment between measurement and predicted values. 
On the basis of these data, applicable to a wide 
variety of sea conditions and ship types, it is 
concluded that the wave-induced motions and 
stresses of ships are distributed according to the 
Rayleigh distribution, for any given set of steady 
conditions of the sea, ship speed and heading. 


ANALYSIS OF LONG-TERM DISTRIBUTION 
FUNCTIONS 


Background 


It is the main purpose of this section to show 
that wave-induced hull stresses and motions of 
ships may be approximated by a logarithmically 
normal distribution if the environmental condi- 
tions under which the data are secured include a 
wide variety of sea conditions, ship’s headings 
and speeds, such as a ship would be likely to 
encounter over a period of weeks, months or years. 
It is a secondary objective to show that the 
heights of ocean waves may be represented by a 
log-normal distribution. 

When this study was initiated several years 
ago, there were few data available on the motions 
and hull-girder bending moments which ships 
experience in service. It was realized that large 
quantities of data would have to be collected in 
order to provide a reasonably firm basis for the 
objective of this program which is the specifica- 
tion and prediction of the service stresses and 
motions of ships. Statistical methods were to be 
utilized in the planning of the tests and the col- 
lection, presentation, and analysis of test results. 
It was felt desirable, early in the planning stage, 
to have a preliminary idea of the general type of 
frequency distributions which might be applicable 
as well as of the length of time over which meas- 


urements would have to be made in order to estab- 
lish the desired distribution patterns. 

Because the stresses and motions of ships are 
caused by the passage of the waves past the ship, 
it was thought desirable to study the extensive 
wave-height observations that had been collected 
over a period of many years both in the Atlantic 
and Pacific Oceans with the expectation that the 
frequency-distribution patterns of ship response 
would, to some degree, be similar to those appli- 
cable to the wave data. Inasmuch as the wave- 
height data comprise visual observations, the dis- 
tributions will be evaluated in a qualitative 
sense only; that is, it will not be possible to assign 
statistical significance levels. The wave observa- 
tions will be analyzed first, followed by the hull- 
girder stress and ship-motion data. Most of the 
analysis will be concentrated on the hull stresses 
because much more stress than ship-motion data 
have been measured to date. 


Heights and Lengths of Ocean Waves 


Fig. 12 shows a chart on which the several 
weather stations, manned by weather ships, are 
marked. For each of these stations the U.S. 
Weather Bureau has collected observations of the 
significant wave height!? and wave period, re- 
ported every 3 hr, over a period of about 6 years 
in accordance with instructions of the U.S. 
Weather Bureau (8). Detailed study of these 
observations indicate that if the data are plotted 
on log-normal probability paper and if a log- 
normal distribution is fitted to the plotted data," 
very good agreement is obtained. Fig. 13 shows 
the distribution for Station C. This curve is 
typical of the distributions obtained from the 
U.S. Weather Bureau data. Station C is the one 
at which the sea trials of the USCGC Casco (1) 
were held. The essential statistical data for sig- 
nificant wave height are given in Table 5. Com- 
parison of the distributions obtained for each 
year with the sum of the observation over a period 
of about 6 years indicates that one year’s observa- 
tion gives a good approximation to the distribu- 
tion obtained over the longer period of time. 
The details of the comparison will not be given 
here. 

Recently J. Darbyshire (9) has obtained meas- 
urements of wave heights by means of a wave 
meter installed on a weather ship. The meastire- 
ments were made over a period of about a year— 
February 1953 to January 1954—at the North 


10 The significant wave height is the mean height of the upper 
third highest waves. 


11 The fitted log-normal distribution corresponds to the value of 
the mean and variance computed directly from the wave data. 
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The ship experienced about 460 variations per hour. 
The test sample consisted of 230 stress variations. 


Stress in Main Deck Amidship, Kips per square inch 
(1 Kip = 1000 pounds) 


Fic. 10 DisTRIBUTION OF VARIATION IN THE LONGITUDINAL STRESS, Main DeEcx, 
AMIDSHIP, AIRCRAFT CARRIER 


Atlantic weather stations I and J, see Fig. 12. 
Darbyshire reported the maximum wave height 
measured each time observations were made at 
3-hr intervals, while the ship was at sea. The 
visual observations made by weather observers, 
on the other hand, are reported as the significant 
wave height. It will be of interest to compare 
the visual observations with the measurements 
that have been obtained with the wave meter. 
If the hypothesis is accepted that the short-term 
distribution of wave height follows the Rayleigh 
distribution, then the maximum wave height and 
the significant wave height are related by a con- 
stant factor. Thus the distributions of maxi- 
mum and significant wave heights should both be 
of the same type, log normal in this case, and 
differ only in their mean values. The US. 
Weather Bureau data indicate that the standard 
deviation”? of the Log, (significant wave height) 
is 0.62 at Station J and 0.61 at Station I as 
compared to a value of 0.57 for the Log, (maxi- 
mum wave height) for the measurements at 
Stations I and J reported by Darbyshire, Fig. 15. 
The wave-meter data have been fitted with a 
log-normal distribution on the assumption that 
the distribution of the maximum wave heights is 
log normal, Fig. 15. The experimental data 
indicate excellent agreement with the fitted 


12 The numerical values given here refer to wave heights measured 
in feet. 


distribution, well within the accuracy of the wave 
measurements. The latter fact, together with the 
good agreement between the standard deviations 
of the distribution of significant and maximum 
wave heights, support the hypothesis that the 
distribution of wave heights may be approxi- 
mated by Rayleigh and log-normal distributions 
for the short and long term respectively. It is 
realized that the visual estimates of significant 
wave heights obtained by the U.S. Weather 
Bureau may not be accurate estimates of the 
mean value of the third highest waves. But as 
long as these estimates are proportional to the 
actual significant wave height, the form of the 
cumulative distribution (log-normal type) and 
its variance remain unchanged; only the mean 
value would change since log (constant) x is equal 
to log (constant) plus log x. 

Fig. 14 shows the distribution of the predomi- 
nant wave periods observed at Ocean Station B. 
Here again the log-normal distribution appears 
to provide a reasonably good fit. If these wave 
periods are converted to wave lengths by means 
of the formula,!? L = 5.12 T?, where L is the 
length of the wave in feet and 7 is its period in 
seconds, which is applicable to gravity waves 
in deep water, it may be seen that the wave 
length will follow the log-normal distribution in- 


13 The numerical value of the constant (5.12) does not affect the 
type of distribution. 
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CUMULATIVE DISTRIBUTION OF VARIATION IN STRESS FOR Main Deck oF AIRCRAFT CARRIER AMIDSHIP 
(See Fig. 10 for corresponding histogram.) 
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14 error of 2 percent in the process of 
mH Numerical integration. 


Line ® represents the log-normal distribution fitted to the 


The dashed lines denote the limits the | 
plotted values may take on due to an 


observed values of significant wave height which 
were obtained at intervals of 3 hours over a period 
of 4% years. The position of the line was computed 
from the observed data. 


represents the distribution of individual wave heights 
derived from the significant wave heights represented [Ii] 
by line A. ins 
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Fic. 13 DIstRIBUTION OF THE HEIGHTS OF OcEAN WAVES AT WEATHER STATION C, 52° N 37° W, NortH ATLANTIC 
OcEAN 
(This distribution is based on 12,365 observations made over a period of 41/2 yr by U.S. Weather Bureau personnel.) 


volving only a change in scale from the distribu- 
tion of periods, Fig. 14. 

Observations of significant wave heights in 
the Pacific Ocean made by merchant vessels 
have also been studied (13). Although these 
data are not as extensive or systematic as the 
Weather Bureau data for the North Atlantic 
Ocean, the indication is that the cumulative 
distribution functions of significant wave heights 
in the Pacific Ocean are approximately log 
normal. 


The log-normal distribution of ‘“‘significant’’ 


wave heights may be thought of as a distribution 
of E where each value of — corresponds to a defi- 
nite sea state which may be alternatively speci- 
fied by a significant wave height, if it may be 
assumed that the Rayleigh distribution holds for 
the distribution of individual wave heights in a 
given sea state. This assumption appears to be 
justified (11). On the basis of this assumption it 
is possible to arrive at the long-term distribution 
of wave heights (as apart from significant wave 
heights) in accordance with the following ap- 
proach. 
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TABLE 6 SS Esso ASHEVILLE—Basic Data ON WAVE-INDUCED STRESSES FOR THE PERIOD SEPTEM- 
BER 1, 1953 To Apri 3, 1954 


Number of Stress 


Type of 
Measurement 
Actual or 
anded Counts 
Load or Ballast 
Days at Sea 
Subject to Study 
Days at Sea during Which 
Suitable Data were Measured 


psi 


Ex 
Subject to Sampling 


Percent of Time at Sea 
Minutes of Recording Time 


4,000-6,000 


Expansion® 
2,000-4 000 


1,690 
23,406 


55 
762 
8 


L 
L 
L 
L 
B 
B 
B 
B 
L 
L 
L 
L 
B 
B 
B 
B 


*The factors 


3 
590 


and kK, are the factors by which the cycles measured on the oscillogram samples are multiplied in order to estimate the number of stress 


Cycles Falling between Stated Limits (for a 30-Day Period at Sea) 


For a Sample of 2 min Taken Every Hour For a Sample of 2 min Taken Every Fourth Hourtt 


Over 12,000 
psi 
Expansion® 
Factork, 
1,000-2,000 
psi 
2,000-3 000 
psi 
3,000-4,000 
psi 
psi 


4,000-5,000 


1,226 
99,171 


758 
61,315 
1,336 
69,873 
06 
64,620 
8 


5 
61,582 


427 
22,332 


8 
111 


7 
785 6 
70,752 


280 
25,172 


125 
11,238 


K 
cycles that would have been obtained if continuous measurements had been taken over aperiod of 30 days at sea. 


+P eak-to-peak variation. 
ttEvery fourth sample recorded was utilized in this analysis. 


Mathematical Synthesis of Cumulative Long- 
Term Distribution of Wave Heights Based on As- 
sumption of a Log-Normal Distribution of Sig- 
nificant Wave Height. Let the instantaneous dis- 
tribution of the wave heights pertaining to a 
particular sea, be defined by the Rayleigh dis- 
tribution, thus 


2269 
= SS e-2/Ei 
p(x) 5,0 


where!4 E, is the mean square value of all the 
individual wave heights comprising the ith sea. 
Let E = cE” be the significant wave height 


where c is a constant. The distribution of F is 
log normal, thus 


— (log E —u)? 
202 


p (log E) = —F>= 
oO Z 


where o? is the variance, and u is the mean value 
oflog H. Therefore 


4 E; is also equal to the total area under the power spectrum. See 
references (14, 15) for a discussion of the power-spectrum concept. 


P(E) = % bllog E) 


and the probability P that x take on a value < x, 
is 


P 1 —(log E — u)? 
= = — 20? 
(#1) a cE V/ 20 
E=0 
xX 
YPig See Z 
[ yma” BE? dx dE 
x=0 


The second integral is integrable, thus 
E =o 


Hl 


— (log E —u)? 


| 


The integral may be evaluated numerically by 
a summation, thus 


P(x;) = 
0 

— c2z42 2 
l—e £& | dE 


J=N 
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Fic. 14 DisTrRIBuTION OF SIGNIFICANT WAVE PERIODS AND LENGTHS FOR PERIOD, JANUARY, 
1949, To JUNE, 1955, BASED ON 16,060 OBSERVATIONS 


Where f; is the fraction of all observed sea states distribution corresponding to each value of E, 
defined by &;. which defines an observed significant wave height 

The summation also could be carried out by and then summing over-all values of observed 
computation of the individual probability density significant wave heights. Both processes were 
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TABLE 7 SS Esso ASHEVILLE—DISTRIBUTION OF WAVE-INDUCED 
STRESS VARIATIONS FOR A PERIOD OF 30 Days AT SEA BETWEEN SEP- 
TEMBER 1, 1953 AND APRIL 3, 1954 


(Average period of measured stress variation was about 9.2 sec. 


Most frequently observed pe- 


riod was in range 5 to 7 sec; see Fig. 3.) 


Estimated Number of Cycles of Magnitudet Greater than the Indicated Stress 
(For a period of 30 days at sea) 


Average No. Cycles 
P and S 
Load and Ballast 


181, 128 
97,640 
33,292 
12,170 

4,174 
1,139 
566 


188,723 
89,552 
28,237 
11,574 

3,889 


Ship Loaded Ship in Ballast 
Starboard Gage Port Gage Starboard Gage Port Gage 


188,007 
97,667 
33,047 
12,727 

5,679 
1,020 
928 


174,072 
100,381 
30,508 
8,176 
2,161 


173,710 
102,958 
41,376 
16,204 
4,966 
2,239 
621 


“This estimate is based on oscillograms taken for a duration of 2 min at intervals of 1 hr. 


>This estimate is similar to the preceding one except that, for stress values less than 4,000 psi, sample 
oscillograms of 2-min duration, taken at intervals of 4 hr, were utilized in Figure 2. 


{Peak-to-peak variation. 


carried out for the wave-height data observed at 
ocean station C. The method of summation given 
by the last expression is by far the more con- 
venient of the two choices. The computed 
cumulative probabilities for the wave heights are 
plotted in Fig. 13. Inspection of this figure 
shows that the cumulative distribution of wave 
heights derived from the distribution of significant 
wave heights also may be represented by a log- 
normal distribution. 

It is concluded on the basis of this analysis of 
the extensive available wave observations, that 
the distribution functions of both the significant 
wave heights and the wave heights proper of 
ocean waves are approximated by the log-normal 
distribution. This conclusion has not been evalu- 
ated by use of standard statistical tests of signifi- 
cance because the data consist principally of 
visual observations and the accuracy of these 
observations has not been determined thus far. 


Stress Variations 


The stress and ship-motion data given in this 
paper furnish a “‘truncated’”’ or “censored’’ dis- 
tribution!’ thus named because only values larger 
than a specific lower limit are given. Fisher (16) 
has devised a method for determining the param- 
eters o and £ which define the normal (or log- 
normal) distribution from the measured truncated 


15 Truncated at the lower end of the distribution. 
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data. If the analytically expressed distribution, 
thus determined, fits the measured data reason- 
ably well, it is probable that the data follow the 
normal distribution. The method devised by 
Fisher will be applied here, except for the cases of 
the Gopher Mariner and Ocean Vulcan for which 
the degree of truncation is so large that a better 
estimate probably can be made on the basis of the 
estimated average periodicity of stress variation. 
The latter two cases are given only to help sup- 
port the conclusions based on other more exten- 
sive tests. 


Esso Asheville 


The period covered by the data is from Sep- 
tember 1, 1953 to April 3, 1954. Oscillograms of 
strain, taken at hourly intervals for a duration 
of 2 min each, were analyzed in terms of stress 
variations (see Tables 6 and 7) and classified 
according to the magnitude of these variations. 
The measured number of variations were then 
multiplied by a factor to give the estimated num- 
ber of variations which would have been re- 
corded if continuous measurements had been 
made. Application of the Fisher method to the 
truncated data yields the log-normal distribution 
patterns shown in Fig. 16. The plotted points 
in Fig. 16 were obtained by use of the theoretical 
degree of truncation; their scatter from the 
straight lines should give an unbiased measure of 
the applicability of the assumed normal distri- 


Station | 


‘| x Observations of Significant Wave Height 
dirom U.S. Wealden Bureau Dlata 1947- 53) 


Stations | and J 
| Measurements of Maximum Wave 
eight obtained by Darbyshire 


H 
(R ee 9) during period Feb! 
1953 - Jan 1954 


Note: The straight lines are the log-normal 
cumulative distributions computed ; 
directly from the experimental data. ( p> 


Probability of Exceeding Wave Height (Percent) 
Sane of Not Exceeding Wave Height (Percent) 


49.99 


5 6 ¢ BY 10 40 50 60 70 80 90 
Wave Height - feet 


Fic. 15 COMPARISON OF WAVE-HEIGHT DISTRIBUTIONS DERIVED FROM VISUAL OBSERVATIONS AND FROM MEASURE- 
MENTS OF WAVE HEIGHTS AT ATLANTIC OCEAN STATIONS I AND J 


bution, within the range throughout which Let the hypothesis be examined that, for this 
stress data were obtained. vessel and service, the logarithm of the sea-induced 
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TABLE 8 SS Esso ASHEVILLE—COMPUTATION OF CONFIDENCE LIMITS FOR DISTRIBUTION OF STRESS 
VARIATIONS, PORT SIDE, SHOWN IN Fic. 16 


1 Confidence Limits 
Log,, Stress hi q= 1 ( P-0.67 
z) faded 2 
z,* (1-p) x 10 
‘b Tower 
h psi psi 


20,000 -1.3255 1.757 1.980 0.99999 | 0.247 4.5480 4.0540 35,300 11,300 
x 10+ (0.494)t (4.7950) (3.8970) (62,400) | (6,400) 
16,000 -1.2286 1.5095 6.798 0.99995 | 0.153 4.3571 4.0511 22,800 11,100 
x 104 (0.306) (4.5101) (3.8981) (32,400) | (7,900) 
~1.1706 1.3703 1.361 0.009 0.99991} 0.103 4.2491 4.0431 17,900 10,900 
x 103 (0.206) (4.3521) (3.9401) (22,500) | (8,700) 
10,000 4.0000 6.718 0.9995 0.0491 4.0491 3.9509 11,300 8,900 
x 103 (0.0982) (4.0982) (3.9018) (12,700) | (8,000) 


8,000 0.9276 0.8604 1.729 0.150 0.9985 0.0335 3.9366 3.8696 8,600 7,400 
(0.0670) (3.9701) (3.8361) (9,300) | (€,900) 


4,000 3.6021 0.6266 0.3926 | 0.1781 2.3 0.977 0.0126 3.6147 3.5895 4,100 3,900 
(0.0252) (3.6273) (3.5769) (4,200) | (3,800) 

3.0000 0.0245 0.00060 47.0 0.0060 3.0060 2.949 1,010 990 

(0.0120) (3.0120) (2.9980) (1,040) (970) 


*The quantiles of zy are distributed normally with the mean value zy and standard deviation i 1 ) 7 where p is the theoretical value of the 
Zz 
L 
cumulative distribution function corresponding to zy and f(z;) is the value of the logarithmic distribution function corresponding to zy. 
N = 4,597 observations. 


**The distribution function of z, (log of stress on port side of deck) has a mean value ze, = 2.9755 and a standard deviation o, = 0.3167, that is 
1 (z, -2,)* e 
f= en ee f = 0.1004 - 2.1650 (z, - 2,)?. 
or ln 20? 10 ame 


{The figures in parentheses are computed on the basis of the assumption that every fourth observation was an independent random observation, 


Net (4,597). 


stresses is normally distributed. Statistical tests quantiles of samples taken from the hypothetical 
of significance assume that the data consist of population. The latter method has been utilized 
random independent observations. Consecutive to fit confidence limits to the distribution of stress 
cycles of stress and motion are not independent.!® variations measured on the port side amidships, 
Thus it will not be possible to give accurate prob- see Fig. 16 and Table 8. It is evident from an 
abilities for the hypothesis considered. However inspection of Fig. 16 that the experimental data 
if it is assumed, falsely, that all data are random are approximated by the log-normal distribution 
and independent then a better fit will be required with a good degree of probability, especially 
between the assumed log-normal distribution when the accuracy with which the stress meas- 
and the experimental data than is necessary at urements were obtained is considered. To avoid 
the given confidence level. confusion, confidence bands have not been fitted 
One rather simple test of significance is Kol- to those straight lines in Fig. 16 which represent 
mogorov’s statistic (17) which permits one to the cumulative distribution function of the star- 
assert with a specified level of probability that board and of the. average stresses. However, 
the true continuous probability distribution func- visual inspection indicates that here too there is 
tion is contained within the confidence band good agreement. 
everywhere. A more accurate (and tedious) test Finally, the hypothesis may be tested by in- 
is to compute the expected variation of the quiring into the relationship that exists between 
the underlying cumulative distribution of all 
16 Autocorrelations have been computed for typical stress and Stress variations and the cumulative ‘‘extreme 


pitch variations; these computations indicate that observations sepa- ” : : : Oa si 
rated by about 6 cycles are ‘‘reasonably”’ independent. value’”’ distribution of the stress variations. 
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TABLE 9 HUuLL-GIRDER STRAIN VARIATIONS MEASURED ON STRENGTH DECK OR KEEL AMIDSHIPS 


cs © 
_ = 9a at 
2/23 2g 3 Be | Sac Frequency Distribution of Stress sos 
a oo = = os ss > 2 
3 [2e]s é SE | oSen — 2 Sp 
s2 > GI = os a.) 2 a a iS 
iS ao S 3 3 3 ta 2ssze Number of Stress Variations per 30 Day Period of Sea Duty 8 g Z 
7 =o i SS SS 
ce rs TH 3 = = >i (The class limits are given in kips per square inch) ms 
ula F>o1 | 
-— =F T 
Emrd 
es Yel 1-2 23 3-4 4-5 5-6 6-8 8-10 | 10-12] Class 
wy ro) Limits 
a] De} 144! 1sep53-3Aprsa] - 220,000 lama 
o>|oc 83,488 64,348 | 21,122 7996 3035 573 505 48 13 
BE) 6 18. 
2 re = Average of Load and Ballast Conditions This distribution is based on 10,580 measured stress variations. The average of the port and starboard 
a gage data was used. 
= —— 2 ag Class 
v7 = : 8-10 10-14 | >14 | Limi 
2 3 é 74 Dec 52 - Apr 53 x wh ERIS 
2 Be Oct 53 - Nov 54 121,706 39,215 2381 176 39 13.9 
| es BBO SE CH 
a ; 
g a= M This distribution is based on 403,343 moasured stress variations 
a Pe | 
paler? 2-4 4-6 8-10 10-12 | 12:20 | >20 | Class 
=| ee Limits 
= oo 
az} Cg | 50/5Feb54-15May54] 7Sec | 362,000 z 7183 100 19.5 1.8 | 0.65 97.6 
os a S This distribution is based on :7,747 measured siress variations. The average of the port and starboard 
- gage data was used. | 
| | alles i Lal 
za| = 2.2-4.5 | 4.5-6.7 | 6.7-9.0 | 9.0-11.2] 11.2-13.4] 13.4-15.7] >15.7 aoe 
SSI eo - 2% Years 8Sec. | 324,000 LES S426 
o>] > 16,690 | 773 «| 46 3 0.17 0.08 0 
OoF>la =| 
This distribution is based on 210,177 measured stress variations. 
a £ T 
= | 5 2-4 4-6 6-8 B10 | 10-20 | >20 Class 
= |2s | 46] Winter 55-56 = = | Limits] ¢4 
[R = 3 109,260) 40,660 2000 1213 127 0.3 
S |} a= This distribution is based on 235,010 measured stress variations. 
rar aot 
eGl| se Yel 1-1 14-2 2-2% | 2-5 55 Class 
=| =s | 17] Winter 55-56 = - Limits) 16.6 
eeles 62,122 | 31,733 13,126 5868 3198 65 
IS Salle This distribution is based on 65,797 measured stress variations. 
= T 
= 0.67.34) 1.34-2.01 ae 2.68-3.35 | 3.35-6.7 | >6.7 class 
ot — 
= 5 61 Winter 54-55 = = 64,603 34,617 | 9444 5365 14.8 
Sz 1 
S35 r This distribution is based on 236,905 measurements of stress variation. | 
JL 


* The percentage of truncation was computed by the method proposed by R.A. Fisher (see chapter 6 of Reference 12) for all cases except the GOPHER MARINER and the OCEAN 
VULCAN. For the latter ships the percentage of truncation was estimated on the basis of an average period of stress variation. 


It has been shown in the literature (6, 19), 
that, if a distribution is normal or log normal, 
the largest values in repeated large samples from 
this distribution have a cumulative distribution 
F(x) of their own (18, 19, 20) which approaches, 
as the sample size becomes larger and larger, the 
form 


F(x) = exp [—e-@-™)¢] or 
F(y) 


where y = (x — m)q@ and is called the reduced 
variate. 

In the present case x is a possible value of the 
extreme of the stress, and ~ and a are constants 
for a particular distribution. Special co-ordinate 
paper has been devised, Fig. 17, on which such 
a distribution will plot as a straight line. Gumbel 
(19) has devised a method for estimating the 
parameters u and a of the distribution of ex- 


exp [—e-]..[1] 


tremes on the basis of measurements of extremes. 
The method is applied here to the extreme sea- 
induced stress variations measured in 2-min 
samples taken at hourly intervals during 19 
trips, each averaging 121 hr underway. 

The resulting extreme-value distributions for 
the port and starboard strain gages are plotted in 
Fig. 17. If all values of x are different, F(x) for 
each observed point equals z/(m + 1) where 2 is. 
the order number of the observed value when they 
are arranged in order of increasing magnitude and 
n is the total number of observed extremes. The 
plotted points represent the individual meas- 
ured extremes whereas the straight line in Fig. 17 
represents the best fit to the measured data, 
utilizing the method of Gumbel (19) under the 
assumption that the underlying distribution is 
of the log-normal type. 

On this assumption, F(x) represents the frac- 
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The straight line represents the 
log-normal distribution computed 
directly from the measured data. 


SS GOPHER MARINER to the Data ——-——-——-—— 


Percent of Variations Larger Than the Indicated Value 


Percent of Variations Less Than the Indicated Value 


lo Pe 


P= 


0.5 0.6 0.7 0.80.9 1.0 
Stress Variation — Kips/Inch2 


0.2 


3.0 4.0 5.0 6.0 7.0 8.09.010 


Fic. 18 CuMULATIVE ‘‘LONG-TERM’’ DISTRIBUTION OF WAVE-INDUCED, HULL-GIRDER STRESSES FOR SEVERAL SHIPS 


tion of all trips which, on the average, are ex- 
pected to give an extreme value of the stress less 
than x. On the bottom scale, values of 1/(1 — 
F), denoted by 7, are shown; T is the number of 
trips the ship must make, on the average, in order 
to exceed the limit during at least one of the 2- 
min sampling periods. The dashed confidence 
limits drawn in Fig. 17 denote the limits within 
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which 67 per cent of the experimental data would 
be expected to fall if measurements of extremes 
were conducted indefinitely and if the assumed 
distribution were valid. Inasmuch as the scatter 
lies well within these limits, it is concluded 
that, for the service conditions encountered, the 
extreme-value distribution (from hourly 2-min 
samples) of the sea-induced hull-girder stresses 


TaBLE 10 SS GopHER MARINER—EXCERPTS OF DATA FOR SEVERAL SEVERE SEA CONDITIONS 


Direction Rel 
to Ship 


Wind and Sea 
Observed Cond. 


Locat 
ey Fully Developed Under 


Computed Data for Waves 


Number of Cyclic Stress Variations Falling Between Given Limits 


= 
= (8 
~~ A e @ |3 eo 
pprox Est'd Observed Conditions 2 Stresses Were Measured at Main Deck Near Amidships 
2 BW o9|-2 
a = SS Sr = aE 
= E/2|s Bl Bl ele gle} 2) 28 | =|sslEe 
SB 5 = = a = 5 3 a) ba = = = a = |a-]ce2.o]4000 psi to | 6000 psi to | 8000 psi to |10,000 psi to| 12,000 psi to Over 
5 2 $ 3 = 2 £| o|Se 2 3 € & 3 = EC\ES! 6000 psi 8000 psi 10,000 psi | 12,000 psi 20,000 psi | 20,000 psi 
2 S |e] Ss BS) SSeS Sa] = | sieees 
g PT alse Sl Seis 8/2] 8) 2] +s [s T 
3 eo/-]5 = QD = |= |Port|stbd | Port| Stbd | Port | Stbd| Port |Stbd | Port |Stbd | Port] Stbd 
Epes) 0600 fe Pott 0-37 | 1400 12.2 3e4| sz | 6 ]4i5 | o |33 | 0 | o | o Jo | of o 
W| 2 | 25-26 | 0400 aa ‘ 1339] 1595 | 109 |181 | 22 |18 | 0 | 3 0 | 0 o|o 
{ 
wi3l 4-5] 1245 = e gig | o |157 8 | 3 | 12 O73 0 o | 
8- 9 | 1130 ‘ 4 1819 | 83 |314 | 21 |35 | 3 | 20 0 | 0 | 0 
W] 4 } 13-14 | 1000 x 7 446 | 8] 6 | PWG WO. 1B OMe ial 0 1 


Notes Taken from Ship’s Log 
12100, 13 Feb to 0100, 14 Feb, Greenwich Time: 
2000 to 0800, 26 Feb, Greenwich Time: 


water over bow and stbd side. 
30900 to 1300, 5 March, Greenwich Time: 


Sea State Code as given by Hydrographic Office Publication No. 606-e published 1950 


*Computed from diagrams | and III of ‘The Forecasting of Sea and Swell Waves’’ Naval Weather Service Memorandum No. 135/45, Admiralty, London, 1945. 


**These are the total number of cyclic variations recorded during the given day. 


Vessel rolling very heavily and yawing widely in very high W-NW sea and very high NW swells. 
Very rough and high W-NW seas. Long, heavy NW swells. Vessel rolling and pitching very heavily, yawing moderately and shipping 


Vessel pitching heavily and pounding occasionally. Shipping water over bow and port side of main deck. Very rough westerly seas and 
short high westerly swells. Reduce speed to ease vessel. 


2h S| 


iS approximately of the form given by Equation 
[1] and that the cumulative distribution functions 
of both the underlying and the extreme-value 
distributions are compatible. 

It is concluded, on the basis of these several 
tests of significance, that the basic hypothesis to 
be tested is ‘“‘probably” true; that is, the loga- 
rithm of the sea-induced stresses is, at least approx- 
imately, normally distributed. 


USS Fessenden 


The cumulative distribution function of stress 
variation was obtained from data measured dur- 
ing the periods December 1952 to April 1953 and 
October 1954 to November 1954. The vessel was 
instrumented with a 10-in. mechanical statistical 
strain gage and counter (2). The measured 
data are given in Table 9 and the corresponding 
distribution pattern is shown in Fig. 18. The 
deviation of the plotted experimental points 
from the hypothetical distribution, which was 
computed from the measured data on the as- 
sumption that the stresses do follow the log-normal 
distribution, was within the accuracy of meas- 
urement. Thus it may be concluded, for the 
USS Fessenden, that the distribution of hull- 
girder stresses or bending moments can be ap- 
proximated by a log-normal distribution. 


SS Ocean Vulcan 


This vessel has been subjected to extensive sea 
tests. The general construction of the vessel 
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as well as the statistical strain gage installed 
thereon have been described elsewhere (21). 
Stress variations which were measured over a 
period of 214 years, have been reported in refer- 
ence (7) from which the values given in Table 9 
have been taken. Owing to the high degree of 
truncation of the data, it was thought best to 
estimate the degree of truncation from the meas- 
ured average period of stress variation. The data 
for this vessel appear to follow the log-normal dis- 
tribution, Fig. 18. 


SS Gopher Mariner 


This vessel was instrumented with two TMB 
mechanical strain-cycle gages and counters (2). 
Measurements were made during six consecutive 
transatlantic crossings, during which the vessel 
encountered several days of very rough weather. 
Excerpts of data are given in Table 10 in order to 
give some idea as to the environmental conditions 
encountered. The degree of truncation of the 
data was estimated from the estimated average 
period of stress variation. Fig. 18 indicates no 
disagreement with the assumption that the log- 
normal distribution applies. 


Destroyer 


The cumulative distribution of hull stress ex- 
perienced amidships in the keel of a destroyer was 
obtained from measurements made during the 
winter season ’55-’56 during operations in the 
North Atlantic Ocean. The strains were meas- 
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Frc. 19 USCGC Casco—SratisticaL DISTRIBUTION PATTERNS OF Sure MOTIONS 
(Data were measured during March and April, 1951, in the North Atlantic.) 
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Average Percentage of Variations Greater than the Indicated Values 


Number of Variations Measured 


Lower 95% Confidence Limit47 A [| 
nl Confidence Limit gy | | | || [| 
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Pitch Angle in degrees, Double Amplitude 


Bow Moving Down 


Level 
1 4 8 12 16 20 
Double Amplitude in degrees 
Pitch Angle 


Fic. 19 (continued ) 
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ured by means of electric wire strain gages. The 
measured data are given in Table 9 and the log- 
normal distribution which was computed for the 
truncated data is shown in Fig. 18. The fractiles 
corresponding to the larger, relatively rarely ex- 
perienced stress variations do not check the fitted 
distribution too well. However, inasmuch as 
the sampling period did not cover an entire 
winter season it can be seen that a single storm 
can account for a disproportionate number of 
the relatively rare, high-stress values. 


Aurcraft Carrier 


The stress data were obtained during the ship’s 
winter season °55-’56 operations in the North 
Atlantic Ocean. The strains were measured by 
means of electric strain gages mounted on the 
keel amidships. The measured data are tabu- 
lated in Table 9 and the corresponding log-normal 
distribution is shown in Fig. 18. The fractiles 
of the experimental distribution agree unusually 
well with the computed log-normal distribution, 
well within the experimental accuracy of the 
measurements. Thus it may be concluded, for 
this ship, that the hull-girder stress or bending- 
moment variations can be approximated by a 
log-normal distribution. 


USCGC Unimak 


The stress data were obtained during the 
vessel’s operation as a weather ship in the winter 
season ’54-’55. The log-normal distribution, 
fitted to the truncated experimental data, is shown 
in Fig. 18. Excellent agreement between the 
plotted experimental fractiles and the computed 
log-normal distribution is indicated. 

It is concluded on the basis of this analysis of 
extensive measurements of hull-girder stresses 
obtained on seven ships of widely different size, 
type, and service that the long-term distribution 
of hull-girder stress may be approximated by a 
log-normal distribution. 


Ship Motions 


The long-term distributions of ship motions 
have been studied for the USCGC Casco and the 
Esso Asheville. 

USCGC Casco. Cumulative distributions of 
heave, pitch, and roll motions taken from refer- 
ence (1), are shown in Fig. 19. It is apparent 
that the logarithmically normal distribution fits 
the data quite well. 

Esso Asheville. Tables 11 and 12 give the 
variations of pitching and heaving accelerations, 
classified according to magnitude and frequency 
of occurrence for both the ballast and load con- 
ditions of the ship. The Fisher method for treat- 


TABLE 11 SS Esso ASHEVILLE—Basic SHip-MOTION DATA FOR THE PERIOD SEPTEMBER 1, 
1953 To APRIL 3, 1954 


Number of Cycles Falling between Stated Limits® (Peak-to-Peak) 
for a Period of 30 Days at Sea 


For Sample of 2 min every Hour For Sample of 2 min every Fourth Hour 


Minutes of 
H 
P 


Actual or 
Expanded Counts 
Load or Ballast 

Days at Sea 
Subject to Study 

Days at Sea during Which 
Suitable Data Were Measured 
Percent of Time at 
Sea Subject to Sampling 
Recording Time 


Expansion®® 
Factor 
Expansion®® 


Heave 
Accel. 
(Amidships) 


> > 10.25 - 0.50 


1,031 
58,798 


Heave 
Accel. 
(Amidships) 


1,208 
49.12 rae 


T 


18,660 | 951 


949 565 
60,337 | 30,550 


*H denotes Heave Acceleration; its intensity is specified in terms of units of gravity (G’s). 
P denotes Pitch Acceleration; its intensity is specified in terms of an equivalent linear acceleration in units of gravity. This value is equal to 
‘angular ey) ft gravity units 
2 100 ft) + 32.2 — = ———_—__ 
( intad/sec2,  J% (100 ft) + 32 sec 100 ft length 


**The expansion factor has been defined in Table 6 


TaBLeE 12 SS Esso ASHEVILLE—DISTRIBUTION OF HEAVE AND PITCH 


ACCELERATIONS FOR A 30-Day PERIOD AT SEA 
(Average estimated period of measured motions is 9.2 sec taken from stress oscillograms.) 


Estimated Number of Cycles of Magnitude ¢ Greater than the Indicated Value 
(For a period of 30 days at sea during the rough weather season) 


Pitch Acceleration 
os s se 


Average 
Las Loaded Ballast 
105,432 


43,552 | 94,185 
8,426 | 20,086] 22,438 | 26,12 19,611 | 33,848 
1,012 4,428 | 2,566 

1,024 


*This estimate is based on oscillograms taken for a duration of 2 min at intervals of 1 hr. 
**This estimate is based on oscillograms taken for a duration of 2 min at intervals of 4 hr (used for plotting). 
fPeak-to-peak variation. 
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ing a truncated distribution was applied to ob- 
tain the results shown in Figs. 20 and 21. The 
classification used for the motion data (acceler- 
ations) was coarser than that used for the stress 
data on this ship. Kolmogorov’s statistic was 
the only measure of significance applied in these 
cases. An inspection of Figs. 4 and 5 indicates 
that the log-normal distribution fits the data 
fairly well. 


APPLICATIONS OF STATISTICAL DISTRIBUTION 
PATTERNS 


A knowledge of the distribution of wave- 
induced ship motions and stresses together with 
the height and length of ocean waves is of con- 
siderable interest in many applications such as 
engineering specification and design and oper- 
ational studies of merchant and combatant ships. 
Following is a discussion of a few representative 
examples illustrating some problems the solution 
of which may benefit from this study. 


Specification of Engineering Requirements 

There are numerous occasions where it is de- 
sired to specify the environmental conditions to 
which shipborne equipment, personnel, and struc- 
tures are subjected. In many cases the ship 
motion will be among the more important items 
that need to be specified. Several methods of 
presenting such information have been given in 
this work; for example histograms, basic and 
cumulative probability distributions. 

In the design of a ship-stabilization system it is 
necessary to estimate the required capacity of 
pumps, motors (22), and so on. The cost of the 
installation is very much a function of the required 
capacities. The designer must therefore be able 
to predict the percentage of time that the roll 
angle will be less than any specified value in order 
to evaluate the desirability of the installation of 
the proposed stabilizer. The distribution pat- 
terns given here are directly applicable to this 
problem. 

In the design of fire-control equipment, air- 
craft launching and landing apparatus, it may be 
necessary to predict the maximum ship motion 
that may be expected under given sea conditions. 
Although it may not be practical to design for 
operation under rare, extremely severe conditions 
it is desirable to have an estimate of the condi- 
tions under which satisfactory operations can be 
expected. The distributions of significant wave 
heights and wave lengths together with the dis- 
tribution of ship motions for various combinations 
of sea condition, heading and ship speed should 
make an evaluation of the operational possibili- 
ties practicable. 
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Often it is desired to estimate the extreme value 
of stress or motion that a ship structure may be 
expected to experience over a specified period of 
time. For example, the extreme value of the 
wave-induced, hull-girder stress that can be ex- 
pected over the service life of a ship is of interest 
inasmuch as the safety of the ship depends upon 
keeping the total hull-girder stress below the ulti- 
mate strength of the structure. 

Therefore, it will be of interest to make an esti- 
mate of the largest stress to be expected in a 
given period of time. As an example such an 
estimate will be made on the basis of the under- 
lying log-normal distribution for the Esso A she- 
ville shown in Fig. 18. For this purpose a formula 
will be developed which will give the fraction f of 
all samples of size NV (belonging to a given under- 
lying distribution) which will have at least one 
(extreme) value of x > xm. The formula is as 
follows 


f = 1 — e-NIl — P(xm)) 


This formula may be developed as follows (as 
was pointed out to the author by Dr. Lieblein of 
the National Bureau of Standards): 

Let P(%m:) be the fraction of the underlying 
(basic) distribution of x with x < xm. 

Take a sample of N values of x in such a way 
that each x is chosen at random from the basic 
distribution, and repeat such sampling many 
times. Then the probability of choosing a value 
less than x,,; when one single value of x is chosen 
from one sample of size N, is P(xmi) and con- 
sequently the probability that some members of 
the sample of size N exceed x 1 is 


to UP Ge) I 


Therefore the fraction f of all (a very large num- 
ber) samples of size NV which will have one or more 
values of x greater than xm, iS 


jal = LP Cally 


This formula is exact, but for large WV it is in- 
convenient. Now for large N it is practically 
certain that in any given sample, Xm, has such a 
value that P(xm,) is close to 1. For values of 
Xm, such that 1 — P(xm,) is small, a more con- 
venient approximate formula can be found. We 
have 


log.(1 — f) = N log. P(%m) = 
N log. {1 — [1 ~ P@n:)}" 
log. (1 — f) w —N[1 — P(%m)] 
by use of the series 


log. 1 = *%) = =x — x7/2.... 


Percent of Variations Greater Thaf the Value Indicated 


Confidence Band® Based 
Motion No. of Measurements on Method on Kolmogorov’s Statistic 


and Load Condition | Which c.d.f. is Based lat a Probability Level of 0.95 
€,,, 0.95=21.3591 


+ 3.30 percent 
+ 3.96 percent 


*The unknown distribution (for the load condition) Is na to lie within the plotted confidence 
band if it is assumed that the measured data were random and ence samples. 
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Fic. 20 SS Esso ASHEVILLE—STATISTICAL CUMULATIVE DISTRIBUTION OF HEAVE ACCELERATION 
(These distributions are based on data given in Tables 11 and 12.) 
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Percent of Variations Less Than the Value Indicated 


Percent of Variations Greater Than the Value Indicated 


Motion 


No. of Measurements on of 
and Load Condition 


Which c.d.f. is Based Sampling 


Pitch Acc., Ballast 
Pitch Acc., Load 


independent samples. 
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*The unknown distribution (for the load condition) is expected to lie within the plotted confidence 
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Fic. 21 SS Esso ASHEVILLE—STATISTICAL CUMULATIVE DISTRIBUTION OF PITCH ACCELERATION 


(These distributions are based on data given in Tables 11 and 12.) 
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Percent of Variations Less Than the Value Indicated 


Miner's Rule for Estimation of ‘Cumulative Damage’’ 0 


ny tae) 


D= + coon 
Ny Ny Ng 


; D=1 Corresponds to failure 


where n; denotes the number of stress cycles applied at stress i 
Nj denotes the number of stress cycles of magnitude i, which if applied alone would cause a fatigue failure. 


S-n C 
Plate-Stiffener 
Combination 


= oo 


Z; 
Y 


Fluctuating Stress as Percentage 
of Static Ultimate Stress 


= 
ss __ 


102 108 108 


n= Number of Stress Cycles per Year (or Load Cyc 


$= 


urve for Typical 


les) 


| 
Total bamage = 18.6% 
Yea 


Life Expectancy = 5.4 Years 


| 


107 
Percent Damage Per Year 


Fic. 22 ILLUSTRATION OF THE APPLICATION OF A LOAD SPECTRUM TO THE DESIGN FOR ENDURANCE 
STRENGTH 


Thus f is given approximately by the expression 
shown in Equation [2]. This result is inde- 
pendent of the form of P. As a first example of 
the application of Equation [2] let it be required 
to estimate the extreme stress variation that 
would be exceeded by one half of all ships” 
(f = ¥%) if 10,400 random stress measurements 
(NV = 10,400) were made on each of many ships. 
Assume that the underlying distribution appli- 
cable to these ships is that given for the port side 
of the main deck amidships in Fig. 16. 


f=%, N = 10400 


From Equation [2] 
f = % = 1 — ¢710400[1— P(zm)] 
—10400[1 — P(*m)] = —0.6931 


or 
1 — P(x) = 6.66 X 10-5 


The corresponding value of the stress variation 
Xm, is read from Fig. 16 as 16,000 psi (from hog to 
sag). Similarly an estimate of the stress vari- 
ation which will be exceeded in 20 winter seasons 
(26,700,000 stress variations) in one half of all 
ships (f = 1) is about 50,000 psi (from hog to 
sag). 

Not much confidence can be placed in the ac- 
curacy of this last prediction because the meas- 
urements on which the prediction is based did not 
include magnitudes of the order of the estimate; 


17 Ships of the Esso Asheville type which would encounter service 
conditions similar to those applicable to the Esso Asheville. 
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i.e., the estimates are not expected to be reliable 
if the estimated extreme value lies far outside 
the range in which experimental data are avail- 
able. 

A better method for making an estimate of the 
probable extreme stress is by use of the experi- 
mentally obtained distribution of extreme values; 
this being the most direct method. 


Prediction of Life Expectancy of Structures 


The results of the study of hull-girder stresses 
(bending moments) given here can contribute to 
the rational design of ship structures on the basis 
of the prediction of the life expectancy of such 
structures. 

The information that is necessary to make such 
a prediction is as follows: 

1 It is necessary to know the service loads or 
stresses that the ship structure is required to with- 
stand throughout its service life. The statistical 
distributions given in this paper provide the wave- 
induced load data in a convenient form. The 
more or less steady stresses such as dead load 
stresses, temperature-induced stresses, and built-in 
stresses also must be estimated. 

2 Adequate data on the endurance strength of 
the ship’s structural components are necessary. 
Probably such data will only be needed for typical 
components as well as for those portions of the 
structure in which a failure. could be disastrous. 

3 Finally, it will be necessary to have a suita- 
ble method or theory available which relates the 
applied service loads to expected structural dam- 


age. Several empirical methods for providing 
such a correlation have been proposed. The best 
known and simplest of these is the cumulative- 
damage law proposed by Miner. 

As an illustration of the manner in which the 
endurance life of a typical structural component, 
say a plate-stiffener combination, can be pre- 
dicted, consider Fig.22. For the sake of simplicity 
it will be assumed that the mean stress level in the 
structure is zero. The histogram of stress varia- 
tions shown in the left of the figure is also known 
as a load spectrum. The damage due to each 
component of the load spectrum is computed and 
the result is plotted as a damage spectrum. For 
example, 100,000 applications of a stress variation 
occur in the range of 4 to 8 per cent of the static 
ultimate stress. According to the S-N curve 
approximately 1 million applications of this stress 
are required to cause failure. Therefore, accord- 
ing to Miner’s criterion, the damage due to the 
application of 100,000 applications of this stress is 
100,000/1,000,000, equal to 10 per cent. The in- 
tegral of the damage spectrum over the stress 
range experienced in service gives the expected 
total damage during 1 year. In the illustrative 
example the expected life of the structure is 5.4 
yr. Owing to the variability of the endurance 
strength of nominally identical structural com- 
ponents, as well as the variability of the service 
loads for structures which presumably are sub- 
jected to the same service, it is necessary to re- 
duce the expected life of the structure by a factor 
in order to arrive at a “‘safe life’’ which may be 
expected to be attained with an acceptably high 
degree of probability. 

The illustrative example was applied to a simple 
case in which the mean stress level was zero. The 
method can be extended easily to the more real- 


istic case where the sea-induced stresses are super-- 


imposed on various mean levels of stress. In the 
latter case it is only necessary to sum the damage 
resulting from the load spectra associated with 
each mean level of stress. 

An interesting exposition of this general ap- 
proach to the design of airplane structures has 
been given by W. Tye (23). 


Prediction of Optimum Operating Conditions 


In ship operations it is often necessary to 
select a heading of the ship relative to the direc- 
tion of the waves and a ship speed such that the 
pitching or rolling motions are below a level be- 
yond which certain operations would become im- 
practical. Examples of such operations are fuel- 
ing at sea and the launching and landing of air- 
craft. For this purpose it is convenient to have a 
graphical or numerical presentation of a parameter 
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Fic. 23. TypicAL PRESENTATION OF THE STaTISTic E 
WuicH DEFINES THE PROBABILITY DISTRIBUTIONS 
(Statistic E is integral of power spectrum. Maximum expected 
magnitude of pitch angle is proportional to VE.) 


which, at a glance, indicates the relative severity 
of motion to be expected for any given ship speed, 
heading, and sea condition. 

Accepting the hypothesis that the short-term 
distributions of ship motions and stresses can be 
represented by the Rayleigh distribution, then 
one can plot the parameter /, which defines this 
distribution, as a function of ship heading, ship 
speed, and sea state. Fig. 23 illustrates such a 
plot. The most probable maximum value of the 
ship motion expected in a given number of oscilla- 
tions is proportional to »/E. Table 4 gives the 
number by which +/E£ is to be multiplied in order 
to give the magnitude of the expected maximum. 

Tables or graphs could be prepared for any 
given ship type to cover the range of operating 
conditions likely to be encountered in service. 

The foregoing examples are illustrations of some 
useful applications of the statistical distributions 
discussed in this paper. 


DISCUSSION 


It has been the main objective in this paper to 
show that the distribution of wave-induced varia- 
tions in the pitch, roll and heave motions of ships 
as well as the wave-induced,'!® hull-girder stress 
variations may be represented by a Rayleigh dis- 
tribution for steady conditions of the sea, ship 
speed, and course, whereas these responses of the 
ship to the sea will be approximated by a loga- 
rithmically normal distribution when distributions 
are considered which extend over a period of weeks, 
months, or longer, during which period the ship is 
subjected to a wide variety of operating con- 
ditions. 

A secondary objective has been to show that the 
foregoing hypothesis is also applicable to the dis- 
tribution patterns of the heights and lengths of 
ocean waves. 

Over 100 individual short-term distributions of 
ship motions and stresses obtained for a wide 
range of environmental conditions on the USCGC 
Unimak were analyzed. A smaller number of 
similar short-term distributions of stresses and 
motions for a large aircraft carrier and for a de- 
stroyer were analyzed in a similar manner. No 
significant difference was found between the hy- 
pothesis and the test data. On the basis of the 
statistical tests of significance, the quantity and 
the scope of the data, one may accept confidently 
the hypothesis that the ‘“‘short-term’’ distribution 
of ship motion and stress approximates to the 
Rayleigh type. 

The long-term distribution of hull-girder stresses 
or motions was determined for a wide variety of 
ship types and operating conditions. The ships 
studied included an oil tanker in coastwise service, 
two dry-cargo vessels, a Coast Guard weathership, 
a destroyer escort, an aircraft carrier, and a de- 
stroyer. Extensive statistical analysis and tests of 
significance were applied to the results of the 
Esso Asheville tests. In this case, it was possible 
to study the extreme-value distribution as well 
as the basic distribution of hull stresses which were 
arrived at by independent methods. The distri- 
bution patterns evidenced by these two distribu- 
tions complement each other; that is, if the basic 
distribution is log-normal then the extreme-value 
distribution of the stress should be of the type 


F(x) 


exp [—e7?] 


The confidence limits applied to the cumulative 
distributions as well as the chi-square test applied 
to the grouped measured data indicate that the 
hypothesis may be accepted with confidence; 


18 It should be noted that the vibratory response of the ship inci- 
dent to slamming is not included in this analysis. 


38 


i.e., the wave-induced stresses and motions, over 
the long term may be represented by a log-normal 
distribution. 

Finally, from a utilitarian point of view, one 
also may accept the hypothesis that the short- 
term distribution of wave heights follows the 
Rayleigh pattern and that the long-term distri- 
bution of wave heights and wave lengths follows 
the log-normal pattern. The latter conclusion 
cannot be evaluated by the use of standard statis- 
tical tests of significance since most wave obser- 
vations were made by eye and no measure of the 
accuracy of these observations has been deter- 
mined thus far. Nevertheless the very extensive 
data obtained by independent observers, the con- 
sistent accommodation of these data to the log- 
normal pattern, and the fact that the response of 
the ship to the waves may be represented by a 
log-normal distribution, all tend to the conclusion 
that the log-normal distribution is applicable to 
the specification and prediction of wave heights 
and wave lengths. 

The existence of a predictable distribution pat- 
tern will greatly simplify the collection and analy- 
sis of service data, as well as furnish a concise and 
graphic manner for the specification and predic- 
tion of ship stresses and motions for given ships 
in a given service. Prediction of the largest 
stresses and motions to be expected over a speci- 
fied period of time can be made readily, although 
the prediction of extreme values is basically less 
reliable than a prediction of those values which 
occur more frequently. The method for predict- 
ing extreme values eventually breaks down by 
predicting values which are too large. This occurs 
because the theoretical distribution cannot be re- 
lied upon at the extreme ranges of the function. 
For example, the theory could predict stresses in 
excess of the ultimate strength of a structure, ob- 
viously an invalid prediction. Rejection of pre- 
dicted values lying outside the range of experience 
would provide a proper safeguard against too ex- 
treme an estimate. 

The parameters computed for the distribution 
functions discussed in this paper are not affected 
appreciably by variations or scatter of those 
values or events which occur relatively infre- 
quently. For this reason, provided the type of 
basic distribution pattern has been established, it 
is possible to obtain a reasonably good approxi- 
mation to the distribution pattern with a limited 
amount of experimental data, because it is not 
necessary to continue measurements until the 
more extreme and infrequently occurring members 
of the population have been sampled. Inasmuch 
as ship motions and hull-girder stresses are of 
rather small magnitudes during the greater por- 


tion of time, it is necessary to measure the smaller 
magnitudes in order to have a sufficiently repre- 
sentative sample for the statistical analysis. Too 
large a degree of truncation of the data would 
make it difficult to establish the distribution pat- 
terns with acceptable accuracy. 

In engineering application very large values 
(extremes) of stress and motion, although they 
occur very seldom, may be of considerable im- 
portance. The reliability of predicting the mag- 
nitude of these extreme, but rare, values is of ne- 
cessity very poor; nevertheless it would be de- 
sirable to have experimental data in this range if 
such predictions are to be made. In any case, 
if statistical predictions are to be made, it is neces- 
sary that measurements of the lower magnitudes, 
which occur relatively more often than the larger 
values be available. However, if rare extremes 
are not of particular interest, the reliability of 
theoretical predictions based on experimental 
data which do not include these rare extremes will 
not be affected appreciably by lack of observations 
of these extreme values. 


SUMMARY 


The analysis of the voluminous data accumu- 
lated during the past several years has indicated 
that the pitch, roll, and heave motions of ships 
as well as the hull-girder stresses, follow the same 
general distribution pattern as do the heights of 
ocean waves. 

It has been shown that the frequency distribu- 
tions of pitch, roll, and heave motion of ships, as 
well as the over-all ship-girder stresses induced by 
the sea, can be represented by the one-parameter 
Rayleigh distribution when the environmental 
conditions of the sea, ship speed, and course are 
constant. Furthermore, the ship responses can 
be represented by the two-parameter logarithmi- 
cally normal distribution when the environmental 
conditions are allowed to vary over the range that 
ships encounter over an extended period of 
service. 

These logarithmically normal distributions thus 
are the result of the physical summation of a large 
number of Rayleigh distributions, each of which is 
weighted in accordance with the relative proba- 
bility with which the environmental conditions 
pertaining to it are encountered in service. 

It also has been indicated that the long-term 
distribution of wave heights and wave lengths 
may be adequately represented by the log-normal 
distribution. A study of wave observations taken 
regularly over a period of 6 years at 14 stations 
in the North Atlantic shows that distribution 
based on the observations taken over 1 year differs 
little from that obtained over 6 years. This fact 
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suggests that inasmuch as ship motions and 
stresses are induced by the waves, 1 year is a 
sufficiently long test period for establishing a 
reasonably good statistical distribution pattern of 
ship motions and stresses. 

The patterns which have been found to fit the 
experimental data can be specified in terms of one 
or two numbers which provides a considerable 
simplification of the problem when it is required to 
state what a ship is expected to do. The presen- 
tation of the Rayleigh distribution has been simpli- 
fied by utilizing a change of variables, such that 
all Rayleigh distributions may be represented by 
means of a single straight line, as in Fig. 3. In 
addition to the basic distribution pattern, the 
form of the distribution of the extreme values of 
the ship response to the sea also has been specified. 

It is desirable, from the standpoint of savings in 
time, and cost, that the actual distribution ap- 
plicable to particular ships be determined in the 
future by means of either theory or model tests. 
There are far too many unknowns to make the 
purely theoretical approach practical. Model 
tests would seem to offer the most direct method 
that could be developed into a useful tool within 
the immediate future. It will first be necessary 
to compare the distribution obtained from model 
tests with those obtained from full-scale sea trials 
under comparable conditions. The required full- 
scale data are at hand and it should not be too 
long before preliminary comparisons can be made. 
The method of St. Denis and Pierson for com- 
puting the power spectrum of the ship response, 
given the power spectrum of the sea as well as the 
speed and course of the ship, is directly applicable 
for the calculation of the frequency distribution 
of the ship motions or bending moments, pro- 
vided certain simplifying assumptions made in the 
derivation of their method are acceptable. Their 
procedure should be of great help in making it 
possible to use simplified model test results in regu- 
lar waves for checking the validity of the model 
test procedure for predicting full scale distribu- 
tions. 

There are many practical applications of the fre- 
quency-distribution patterns. In the determina- 
tion of the capacity of ship and shipboard stabili- 
zation equipment, it is necessary to have a reliable 
estimate of the probability of exceeding given 
angles of roll or pitch. The design of aircraft 
landing gear, rocket launchers, and fire-control 
apparatus requires a knowledge of the ship mo- 
tions expected in service. The ability to land 
planes on a carrier in a given sea can be predicted 
on the basis of frequency distributions such as 
given herein. 

In the design of the ship structure, endurance 


strength may well be a factor. Alternating 
stresses of moderate magnitudes superimposed on 
high mean-stress levels do reduce the strength of a 
structure subjected to repeated loading; further- 
more, in places where cold-working is involved as 
a result of such stress combinations (as appears 
likely near welds) the notch sensitivity and transi- 
tion temperature of the affected material are 
generally raised. In the design of structures for 
endurance strength, the frequency distributions 
here should be directly applicable. 
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Discussion 


Mr. E. V. Lewis, Member: The author has 
carried out an ambitious program of assembling 
data on ship motions and stresses at sea and has 
demonstrated that statistical methods can bring 
them into some kind of order. It is particularly 
encouraging that he finds the short-term ampli- 
tudes (from peak to hollow) follow the Rayleigh 
distribution. Other writers have indicated that 
this should be so, but the author has made a dis- 
tinct contribution in assembling such an impres- 
sive weight of evidence to prove that the Ray- 
leigh distribution does apply to waves, ship 
motions, and hull stresses—over a wide range of 
ship types and conditions. Thus the great value 
of this paper is believed to be the way in which it 
adds to and fills out a growing knowledge of the 
sea and ship behavior. 

It may help to see better how his work fits into 
the broad picture, by noting, for example, that 
the Rayleigh distribution is the same as that 
referred to as the ‘“‘target or chi-square” distribu- 
tion by St. Denis-Pierson (reference 15 of the 
paper). The author’s £ is in principle the same 
as the R of St. Denis-Pierson—both are pro- 
portional to the mean square value of the record 
or the area under an energy spectrum. But the 
author’s — is numerically 4 times R. It is the 
writer's understanding that the ‘variation in 
pitch angle’ and other quantities plotted—as in 
Fig. 2 of the paper—refer to the successive double 
amplitudes or peak-to-trough distances on the 
record. It should be noted that the Rayleigh 
distribution of amplitudes is entirely consistent 
with a normal or Gaussian distribution of points 
on the record, chosen either at random or at 
equal intervals of time. Since theoretically the 
Rayleigh distribution applies only for “narrow” 
spectra, it may be concluded from the author’s 
work that most sea and ship response spectra are 
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narrow, in this sense. It is only necessary then 
to know one parameter, FE (or for practical 
purposes o”), to define the distribution. 

Model research at the Stevens E.T.T. in irregu- 
lar head seas has shown that E for motions and 
for external bending moment (hence stress) can 
be predicted accurately from a knowledge of the 
sea spectrum and of the characteristic response of 
the model to different wave lengths (frequencies), 
which easily can be obtained experimentally. 
Hence, it is comparatively simple to predict the 
distribution patterns, as noted by the author. 
(Model results on motions were presented before 
the Society last year, and work on bending 
moments sponsored by the S-3 Panel of the Hull 
Structure Committee is to be published soon.) 

Actual distribution curves for amplitudes of 
model response have not been published, and 
therefore it may be of interest to show some 
typical plots we have just worked out for com- 
parison with the author’s full-scale data. 

Fig. 24 shows distribution of wave, pitch and 
heave amplitudes for the Series 60 model re- 
ported on last year; speed 12 knots in high 
irregular sea (5 runs). The wave spectrum and 
the motion spectra were quite narrow and very 
good agreement is shown. It should be added 
that a typical case of 18 knots in a more moderate 
sea—with a broader spectrum which according to 
Pierson and Neumann is not so realistic—did not 
show such good agreement. 

Fig. 25 shows distribution of bending moment 
for T-2 Tanker model reported in 1954. Again 
good agreement is shown at a speed of 1414 knots 
in a high irregular sea (3 runs). 

The author’s finding that long-term sea and ship 
data follow a log-normal distribution is also of 
interest, since it makes it possible to characterize 
a mass of accumulated data by two parameters— 
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the mean log value and E (or co”). Variations in 
the parameters for different series of data, and 
hence in the lines plotted in Fig. 18 for example, 
are presumably due to such factors as the service 
in which the ship is engaged (i.e., average severity 
of the weather), the average loading of the ship, 
(i.e., always full or sometimes light), longitudinal 
weight distribution, and in the case of hull-girder 
stresses, the section modulus and other structural 
design features of the hull itself, and even the 
particular point at which the stress is measured. 
It does not seem surprising that if sea data exhibit 
a log-normal distribution (Figs. 13-15 of the 
paper) the ship responses also should follow this 
distribution. 

It is somewhat disappointing that Fig. 18 does 
not extend to larger P values and hence larger 
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stresses, since there is some indication of depar- 
tures from the lines at the top of the graph. Table 
9 indicates that a number of stresses have been 
recorded above 10 kips/sq in. 

Coming next to the application of distribution 
patterns of motions, there can be no doubt that 
they will be of great value for the design of ship- 
borne equipment, stabilization systems, fire- 
control equipment, and aircraft launching and 
landing devices. As for hull-girder stresses, the 
author has performed a service in calling attention 
to the possible importance of endurance strength 
to the ship-fracture problem, and the distribution 
patterns are a necessary step in studying the 
problem. He is quite right in giving attention 
also to the extreme values that might be expe- 
rienced over the life of the ship. In the writer’s 
opinion, this is one of the biggest problems facing 
us in structural design at the present time. The 
statistical approach to this problem is interesting, 
but as he points out, is unsatisfactory. He says, 
“a better method....is by the use of the experi- 
mentally obtained distribution of extreme values.”’ 
Presumably this means assembling much more 
statistical data, with emphasis on high values, 
and I am heartily in favor of this. However, I 
object strenuously to the statement, ‘‘Rejection 
of predicted values lying outside the range of 
experience would provide a proper safeguard 
against too extreme an estimate.’ This simply 
begs the question, and if such a policy were fol- 
lowed inevitably would doom many ships to the 
bottom. 

Even with a great deal of data on high stresses 
recorded at sea, it appears doubtful that a satis- 
factory answer to the problem of expected maxi- 
mum values can be obtained by statistical methods 
alone. There are definite physical limits to the 
bending moment exerted on the hull by the sea; 
wave steepness has an absolute maximum and the 


depth of the ship itself restricts the severity of the 
buoyancy distribution. But these physical limits 
on the external loads actually may provide the 
key to the problem of maximum stresses, by 
permitting the maximum loads to be determined 
on physical grounds. Hence, it is suggested that 
attempts be made to evaluate the maximum loads 
on the basis of purely physical consider- 
ations. 

In this connection, the suggestions of Hazen and 
Nims in 1940! may provide a promising lead. 
They made the following suggestion, assuming a 
satisfactory method for the computing machine 
calculation of bending moments (a joint project 
of the Society’s S-3 and H-7 Panels is making en- 
couraging progress here): ‘‘With a wave of any 
given shape at any given position along the hull 
assume a position and velocity for the ship in the 
heave and the pitch co-ordinates. “These assump- 
tions may be considered sufficient to fix the forces 
acting on the ship since the unbalanced force and 
moment, and consequently the accelerations, are 
fixed thereby. By investigating a series of such 
conditions with the hull having a wide range of 
positions and velocities with respect to the water 
it is possible that an envelope for moment curves 
would be defined which would represent the maxi- 
mum moments that would be encountered by the 
hull for any position and velocity to which it 
might be subjected...” Once having deter- 
mined the most extreme bending moment physi- 
cally possible, statistical methods could be used to 
determine the probability of occurrence of this 
extreme condition. 

In conclusion, I feel that the statistical analysis 
of sea and ship data presented here is of great 
value to the profession, and the work should be 
continued actively. I also agree with the author 
that simple model tests can continue to be of as- 
sistance in connection with such studies. 


Mr. J. Darpysuire,” Visitor: The author 
quotes data which I published in 1955 relating to 
wave records taken by the O.W.S. Weather Ex- 
plorer with a ship-borne wave recorder devel- 
oped at the National Institute of Oceanography. 
These records were taken at the positions 
“INDIA” (61°00’N,15°20’W) and JULIETT 
(52°30’N,20°00’W) and covered the period Feb- 
ruary, 1953, to January, 1954. I have since 
checked the logarithmic normal distribution law 
by using all the data available from the ship from 
February, 1953, to December, 1955. The distri- 


19 “‘Calculation of Motions and Stresses of a Pitching and Heaving 
Ship,” by H. L. Hazen and P. T. Nims, Trans. SNAME, vol. 48, 
1940, p. 111. 
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bution law still appears to be satisfactory although 
there is some departure from it at heights below 5 
ft and above 40 ft. The relation 


p = 6.6 exp — (log x/9.0)?/0.126 


was obtained where p is the probability density 
expressed as a percentage and x is the maximum 
wave height. 

A similar relation was found for wave heights 
recorded at Casablanca from 1952-1955 and it 
was 


pb = 18.0 exp — (log x/3.5)?/0.187 


A similar relation also could be fitted to the 
maximum wave heights recorded at Perranporth, 
Cornwall, England, during the year 1946-1947 
and it was 


p = 13.0 exp — (log x/4.0)?/0.199 
This work has been published in greater detail.?! 


Mr. J. P. Comstock, Member: The two objec- 
tives of this paper as stated by the author ob- 
viously have been attained. 

I concur with the author that ‘there are many 
applications in which a knowledge of the fre- 
quency distributions of hull motions, stresses, and 
the heights of ocean waves can be used to advan- 
tage.” However, when it comes to deducing 
“the extreme value of the wave-induced stress that 
can be expected over the service life of the ship,”’ 
I am not yet convinced that a study of statistical 
distribution of stress can do this for us. 

To apply the example of estimating the safe life 
of a structure to a ship’s hull would require the 
tacit assumption that ship failures are basically 
fatigue failures. Since, in all the voluminous rec- 
ords of ship casualties, there has been found, so 
far as I can determine, no correlation whatever be- 
tween age of ship and incidence of fracture, I do 
not believe that ship fractures are basically fatigue 
failures. I think that in most cases they result 
from isolated fortuitous combinations of severe 
conditions which in all probability most ships 
never experience in their entire life, but which 
nevertheless, if they could be known, should be the 
basis of hull structure design. 

The author states, “In engineering applications 
very large values (extremes) of stress and motion, 
although they occur very seldom, may be of con- 
siderable importance.” So far as hull-girder 
stress is concerned, this is a gross understatement; 
they are of controlling importance. Yet, as the 
author warns us, “‘theoretical distribution cannot 


21 ‘Distribution of Wave Heights,” by J. Darbyshire, Dock and 
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be relied upon at the extreme ranges of the func- 
tion.” This forces us back, in merchant ship de- 
sign at least, to the basis of long-range service 
experience. 

This is not to deny or even belittle the other ap- 
plications of this sort of study suggested by the 
author. 


Mr. Joun Vasta, Member: The author has 
demonstrated in a very able and convincing pres- 
entation that ocean waves and ship response, 
both as to motion and hull-girder stresses, fit into 
the framework of a Rayleigh frequency distribu- 
tion. This is a novel and worth-while contribu- 
tion to our knowledge of the behavior of ships at 
sea. One might logically ask now: Of what prac- 
tical significance is this new information to the de- 
signer? With today’s emphasis on special pur- 
pose ships, which impose exacting environmental 
requirements, the determination of the motions 
and structural deformations a ship is likely to ex- 
perience in a particular seaway is of extreme im- 
portance. The question however still remains: 
Can we predict the ship response from known 
theories during the design stage or do we have to 
build the ship first and then send it out at sea to 
observe its behavior? The author states that the 
actual Rayleigh distribution applicable to a par- 
ticular ship can best be obtained by appropriate 
model tests. He does not consider that this can be 
done by theory alone in its present form. Since 
model tests are time-consuming this means that, 
for planning purposes, adequate lead time must 
be allotted in the early design stage to obtain the 
optimum hull shape for ship motion without hav- 
ing compromised the speed characteristics ad- 
versely. 

From the hull-girder strength standpoint, Fig. 
18 is particularly significant and most illuminat- 
ing. It shows for example that for a destroyer 
hull type a 10,000-psi stress range or less occurs 
99.99 per cent of the time that the ship is at sea, 
and conversely, a stress range of this magnitude 
is exceeded but 0.01 per cent of the time. It. is 
realized that this stress range considers sea loads 
alone. However, since the design stress range 
for this type of hull is well over 30,000 psi does it 
follow that the destroyer’s structure reflects a 
disproportionate conservatism, or does it mean 
that the margin between 30,000 psi and 10,000 psi 
stress range constitutes a judicious insurance 
against those rare and extreme sea loads that 
occur but 1/100 per cent of the time? Obviously 
economies in structural design can only follow 
after a realistic evaluation of the sea loads is made. 
Does the author feel that our present hull struc- 
tures are designed too conservatively? 
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An important feature in the paper is Fig. 22 
which illustrates a rule for estimating the ‘“‘cumu- 
lative damage’’ factor for any typical structural 
component once the fluctuating stress pattern for 
the component is known. The life expectancy of 
5.4 years has been worked out for a hypothetical 
case. It would have been much more significant 
if the example could have reflected a realistic case. 
The choice of the numbers for the purpose of this 
illustration is somewhat unfortunate for it creates 
the feeling that the useful life of ships at sea is 
relatively short, whereas the records indicate 
there have been many successful ships that have 
plied the oceans for 25 and 30 years without hav- 
ing had to worry about their endurance limits. 

The author’s contention regarding life expect- 
ancy of a ship with endurance limit as the con- 
trolling factor is misleading. It is general knowl- 
edge that ships have broken in two under still- 
water conditions where the question of fatigue of 
metal was not involved. Moreover it has been 
established that by using a steel with a lower 
transition temperature a major improvement in 
life expectancy results even though endurance 
limit and stress levels are not affected. Thus, al- 
though the author is careful to acknowledge the 
“variability of the endurance strength of nomi- 
nally indentical structures,’ it is the writer’s 
opinion that the variability is of much more prac- 
tical significance than the endurance-limit calcula- 
tion. 

In discussing prediction of life expectancy of 
structures, the author suggests that the quasi- 
steady stresses, such as dead-load stresses, tem- 
perature-induced stresses, and built-in stresses 
must be estimated, and added to the sea-load 
stresses in order to reach a more realistic approach 
to the ship structural design. This suggested 
course should be approached with considerable 
caution. The writer contends that before new 
design procedures are instituted, one must balance 
the additional expenditure of efforts dictated by 
the new approach with the returns that he expects 
to obtain. In the final analysis the most effective 
design procedure is that one which is simple and 
direct to apply, and which in addition has had the 
benefits of accumulated satisfactory sea-service 
experience. 

One can speak with ardor about combining 
several load factors such as thermal stresses, weld- 
ing residual stresses, and so on, to the sea-induced 
stresses. The truth of the matter, however, is 
that if these factors were considered in today’s de- 
sign procedures the resulting hull structure would 
be much too heavy, and most likely economically 
unacceptable. The meticulous designer might 
welcome an elegant rational design approach, but 


the prospective shipowner might frown at the 
rationality of the refinements especially when the 
weight of the hull is increased over and above what 
is safe past practice. The writer wishes to em- 
phasize that improvements in the hull structural 
design are made not so much by acknowledging 
the existence of many complex parameters and 
pleading for their acceptance for design purpose, 
but in determining how important and really sig- 
nificant they might be. Only research can separate 
the wheat from the chaff, and it is here that inves- 
tigations such as represented by the subject paper 
can play their most important role. 


Mr. WILBUR Marks,”* Visitor: In recent years 
there has been a flurry of activity in the field of 
ship motions, hull stresses, and the like. Theoret- 
ical developments, though not fully verified, 
have reached the point where, given a particular 
ship and seaway, certain ship-behavior statistics 
may be predicted. One of the most prominent of 
these theories is mentioned in this paper (15). 
Although simplifying assumptions are used, the 
mathematics take on such a formidable appearance 
as the development progresses that the end result 
of statistically reducing the motion spectra, in- 
stead of providing a handy tool for ship-motions 
prediction, has created rather an atmosphere of 
uncertainty in many places. 

Following this work by St. Denis and Pierson, 
a hue and cry arose, the substance of which was 
the demand for a simplification of the treatment 
and a set of rules (or graphs) from which a ship’s 
captain, for example, could tell immediately his 
vessel’s response for a given set of conditions. At 
this writing no such simplicity has evolved. The 
reason for this will be discussed. The fact re- 
mains that, although the answers are obtainable in 
principle, electronic analyzing equipment or ex- 
pensive calculating machines are required. There 
is little doubt that a prediction system, albeit sta- 
tistical, that is simple to apply is needed. The 
present paper attempts to solve the ship-perform- 
ance problem by utilizing the experimental ap- 
proach. The difficult task undertaken by the 
author can be appreciated readily. If any aspect 
of ship behavior under any particular set of condi- 
tions is considered to be an “‘event,”’ then the in- 
finite population of events being investigated is 
comprised of the different behaviors of all ships, 
in all sea states, traveling at all speeds and head- 
ings. The object of the experimental approach is 
to describe the population by a law (or laws) 
derived from a sample (or samples) whose proper- 
ties are believed to be the same as those of the 
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population. This then is the method used by the 
author in his attempt at solving the problem. 

To be sure, the data collected are voluminous 
and a determined assault on such a veritable moun- 
tain of information is in itself praiseworthy. Yet 
it is implied in the summary that, to answer the 
question for a particular ship under particular 
conditions, more data, either full scale or in the 
model tank must still be collected. The fact is 
that the sample analyzed here is not in itself 
sufficiently adequate to answer all questions, but 
is intended rather to point the way to the solution. 

There is little doubt that the analysis demon- 
strates the applicability of the Rayleigh and log- 
normal laws to the data sampled, even if the possi- 
bility that other laws may fit the data as well, or 
better, is not precluded. The difficulty lies rather 
with the number /. It is this factor which deter- 
mines the slopes of the straight lines which pass 
through the origin and this must be known in 
order to find the distribution of x for any par- 
ticular case. For every possible combination of 
conditions there is an /, and in order to solve the 
seaworthiness problem a representative sample of 
E must be obtained. Just how many are needed? 

Consider for a moment, the number of tests re- 
quired to make up an adequate sample of the 
population. The experimental variables might 
be: (a) 7 ship parameters (hull oscillations, heave 
acceleration, bow submergence, etc.) (b) 7 basic 
ship hull designs, (c) 3 loading conditions, (d) 5 
seaways, (e) 5 ship speeds and, (f) 5 headings. The 
result is that 18,375 separate tests, either full scale 
or in the model tank, must be made, in order to 
characterize the population. Can data collection 
on such a scale be practical? 

Suppose for the sake of discussion such data 
have been collected; what then? It will then be 
required, as it was of St. Denis and Pierson, to 
present these data in a simple form. If graphs of 
ship performance, as a function of ship speed and 
heading, are assembled, there will be 735 curves 
through which a ship’s captain might have to pore 
and even then be faced with the problem of inter- 
polation between different designs, loading condi- 
tions and seaways. It is quite possible that the 
presentation can be simplified and indeed the es- 
tablishment of a central forecast agency, equipped 
with high-speed computers or a punch-card sys- 
tem, may solve the prediction problem but, as 
long as prediction rests with ship personnel, the 
data will necessarily be in a clumsy form, simi- 
lar to those just described. This perhaps is the 
reason why the theoretical solution has not yet 
been reduced to a simple set of rules. By these 
arguments, I feel that the author has assigned 
himself a difficult task indeed, if fie seeks to pur- 


sue this method tu fruition; namely, the predic- 
tion of service conditions of ships. 

I think there is however a very definite place 
for the fine work presented in this paper. It 
seems reasonable that the ultimate solution to the 
ship-behavior prediction problem lies in a com- 
bination of theory and experiment. Theory must 
define the broad general aspects of the problem 
while experiment verifies and modifies the sim- 
plified mathematical model. As an illustration 
of the value of the material here presented, there 
has just been reported by Cartwright and Rydill*® 
the results of an investigation which clearly 
demonstrate the importance of full-scale experi- 
ments. The roll and pitch of a vessel (RRS Dis- 
covery IT) were measured in long-crested seas for 
12 headings and at 614 knots and also head-on 
to short-crested storm seas. Waves were re- 
corded by a ship-borne wave recorder. The rec- 
ords were spectrum analyzed and the amplitude 
distributions predicted. At the same time the 
motions were computed by the simplified equa- 
tions of motion and again by the “‘advanced” 
equations which include coupling terms. Their 
findings are worth repeating: ‘“The practical re- 
sults are compared with theoretical relations 
derived by methods which ignore certain effects 
but are simple to apply, and advanced methods 
which estimate some terms with greater precision 
but at the cost of greater computational labor. 
The theoretical results are in general quite close 
to and never very different from, the practical 
results, the differences largely being due to statis- 
tical variations.” 

To be sure, the results of Cartwright and Rydill 
are encouraging, but nevertheless their experi- 
ments were not comprehensive and certainly not 
conclusive. This paper, however, comprises a 
cross section of material that covers almost every 
aspect of the ship-behavior problem that has been 
investigated theoretically. In addition, the data 
have been reduced to their probability distribution 
functions. If now theoretical predictions were 
made for each pertinent case analyzed in this 
paper and compared with the author’s results, 
the work of Cartwright and Rydill will have been 
carried a good deal further in the direction of veri- 
fication of theory. This appears to be the most 
profitable approach, at this time. 


Pror. B. V. Korvin-KrouKkovsky, Member: 
The author has presented an excellent paper which 
invokes admiration, but also a certain amount of 
criticism. 


23‘‘The Rolling and Pitching of a Ship at Sea,’’ by D. E. Cart- 
wright and L. J. Rydill, Trans. INA, vol. 98, 1956. 
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The positive side is evident: A tremendous 
amount of careful observations were made on 
several ships under many different conditions and 
the results were analyzed-skillfully and presented 
in a compact and clear set of graphs. The very 
complex and varying conditions of sea and ship 
motions are shown to follow a simple statistical 
law on the basis of which certain conclusions and 
predictions can be made. 

On the negative side is the over-concentration 
on pure statistics, and the lack of any reference 
to the physical nature of the situation (except in 
the excellent section on endurance strength). 
This aspect of the problem is taken up in the pre- 
vious literature on this subject which should not 
be ignored. Two papers will.be cited here as in- 
structive examples. 

The first is by A. J. Williams*4 who has shown 
a comparison of the Rayleigh distribution with 
experimental observations on the rolling of a 
ship at sea. The importance of this reference lies 
in that, although the term “Rayleigh distribution” 
was not used, the expression was derived for use 
in this particular case on the assumption of sea 
waves being composed of a number of waves of 
different amplitudes A and phases ¢ traveling in 
one direction. The waves were of a somewhat 
peculiar type in that the amplitude A was as- 
sumed to be a function of time, while the angular 
frequency w was taken as constant. Thus the ex- 
pression was shown to result physically from an 
event which has two degrees of freedom in A 
and @. 

It is clear from the derivation that it applies 
directly to sea waves. In application to a rolling 
ship an additional physical factor enters in, that 
of the ship response to waves. Thus it can be 
seen that in using the Rayleigh distribution 
several factors are neglected, most obvious of 
which are the spectral composition of waves in re- 
gard to the frequency, the well-known variability 
of wave directions at sea and the ship-response 
factor for a specified type of motion (i.e., pitching 
or rolling) and the direction of ship travel with 
respect to the mean direction of wave propagation. 
In view of these physical conditions, it is a pleas- 
ant surprise that Figs. 2, 4, 6, 8, and 10, show 
such a good agreement between the Rayleigh dis- 
tribution and the observed data. Such discrep- 
ancies as are found to exist should not, however, 
be attributed to chance errors, but rather to the 
failure to use a distribution function corresponding 
to the appropriate number of degrees of freedom 
of the wave-ship system under consideration. 

As the second reference on sea waves the paper 


24 “*An Investigation into the Motions of Ships at Sea,’”’ by A. J- 
Williams, Trans. INA, vol. 94, 1952. 
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by Willi Langmaack” can be cited. In the pres- 
ent connection the reference is of particular in- 
terest in demonstrating the value of the ‘‘normal 
distribution” method. As Dr. Jasper points out 
in the early part of his paper, the Rayleigh distri- 
bution is convenient in that it depends on a single 
parameter E. While for many problems the re- 
sultant common form of distribution is of value, 
for others it may be of interest to bring out the dif- 
ferences of existing conditions. In this particular 
case the “normal distribution’ brings out the 
mean value of a quantity, and the spread of dis- 


% ‘‘Haufigkeitsverteilung der Schwingungsperioden und ampli- 
tuden eines Schifies in Seegang,’’ by Willi Langmaack, Werft- 
Reederei-Hafen, Heft 13, 1941, pp. 204-211. 
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tribution about the mean value. Thisis illustrated 
by Fig. 3(0) of the reference mentioned® (Fig. 26 
of this discussion). The plot represents the fre- 
quency distribution of rolling periods plotted 
against the period. Two values of the maximum 
frequency of occurrence can here be seen at the 
periods of 10.7 and 26 sec, with the normal distri- 
bution about each of these. This plot evidently 
represents two superposed oscillations, the first at 
the frequency of wave encounter and the second at 
the natural frequency of ship’s rolling. Aa exam- 
ple of distribution of wave periods is shown on 
Fig. 6(a), * (Fig. 27 of this discussion), where the 
total motion is seen to be the result of three super- 
posed wave trains. 


The two examples taken from the previous liter- 
ature have shown that ship motions in a seaway 
follow certain well-known statistical laws. The 
present author collected much more data and pre- 
sented the case so forcibly that no possible doubt 
about this can exist. On the other hand, it ap- 
pears to be clear that it is not advisable to con- 
centrate completely on one method of statistical 
presentation, as in doing so certain valuable fea- 
tures of wave and ship motions or ship stresses 
can be overlooked. In particular, a distribution 
of a Rayleigh type should not be taken as an arbi- 
trary statistical law of a fixed form, but the ap- 
plicability of a given form to the problem must 
be considered on a physical basis. The distribu- 
tion will take somewhat different forms (i.e., more 
peaked or flatter) depending on the number of 
variables affecting the variate which is under 
study. 


Dr. V. G. SZEBEHELY, Member: The writer will 
begin his comments by quoting Dr. von Karman: 
“The layman is amazed and the expert is left 
wondering.’’ We are amazed by the fact that 
such a great amount of data could be organized 
according to relatively simple basic laws. The 
profession is grateful to the author for collecting 
this useful information and we wish to congratu- 
late him for the clear and orderly presentation of 
the data. The amazement and wonderings of 
the reader might justify the following questions 
in order to evaluate the merits of the paper under 
discussion. 


1 Why log-normal and why Rayleigh distribu- 
tions? Have other laws been tried also? Is there 
a more basic physical law which governs the 
phenomena and from which the foregoing named 
statistical distribution laws can be derived? Does 
the long-term distribution follow from the short- 
term distribution, or to state the question another 
way—should there exist in principle a definite 
relation between the two distributions which, 
after all, represent basically the same phenomenon? 

It is suspected that the answer—to at least 
some of these questions—is affirmative. 


2 How are the results of the paper to be used 
to answer the ever recurring questions: What will 
be the performance of a ship in a given sea? If 
two sets of ship lines are given, which vessel should 
be built for a certain operation? Willa ship with 
a U-bow lose more speed in a state 5 sea than a 
ship with a V-bow? What is the speed at which 
slamming starts for a given design in a given sea? 
What is the threshold for operation; i.e., what is 
for instance the maximum speed and maximum 
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sea state in which the pitch angle will stay below 
5 deg? (It is possible that these questions cannot 
be answered from the impressive amount of data 
in this paper.) 

3 What is the real significance of the line giv- 
ing the hull-girder stress distribution for a destroyer 
in Fig. 9? Would this line be entirely different 
for another destroyer? Do we have to measure 
stresses for long periods of time on several de- 
stroyers in order to find out the relative merits of 
the various designs? Wouldn’t it be simpler to 
run several destroyers parallel in the same sea for, 
say, 30 min and compare the results? 

4 An elaboration on the claim of the paper— 
“to predict service conditions for ships’ —would 
be appreciated. What is the relative merit of the 
following approach? Determine the transfer func- 
tions of a given vessel in regular waves on model 
scale and then, using the energy spectrum of a 
certain sea, find the expected behavior of the 
vessel? 

Answers to these questions will not destroy 
the appreciative amazement and might help to 
guide the wonderings. 


Mrs. MarGaret D. BLEDSOE,” Visitor: This 
writer recently had the opportunity to perform 
statistical analyses of motion records obtained 
during full-scale trials of Dutch destroyers. It is 
believed that these results might serve to cor- 
roborate the work of the paper under discussion. 
The analyses were performed both by using the 
author’s approach of assuming the motions to be 
described by the Rayleigh distribution and by us- 
ing the autocorrelation technique to find the asso- 
ciated power spectrum. The two methods should 
yield the same results for the total energy provided 
most of the energy is concentrated in a narrow fre- 
quency band, i.e., the assumption of the Rayleigh 
distribution is valid. From a knowledge of the 
total energy one can, of course, determine such 
quantities as significant amplitude of motion, 
average height of the upper tenth highest waves, 
and so on. 

The analyses were performed for several samples 
of pitch data where a continuous recording was 
obtained for 10-min intervals. Of the samples 
analyzed, some were obtained in head seas, both 
moderate and heavy, and some were obtained 
in moderate quartering head seas. Ship speeds 
varied from 12 to 22 knots. The theoretical 
Rayleigh distribution was computed and com- 
pared with the experimental histogram and in all 
cases the Rayleigh distribution was found to de- 
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scribe the actual data with a degree of significance 
comparable to that shown by the author. 

By reanalyzing the data according to the power- 
spectrum concept, the total area under the spec- 
trum was obtained and this was compared with 
the author’s parameter E which defines the 
Rayleigh distribution. The maximum discrepancy 
found was of the order of 7 per cent. This, of 
course, means only a 3.5 per cent error in the de- 
termination of the significant or average height 
of the motions. 

One sample might be cited in particular for 
which the test conditions corresponded to those in 
Fig. 6 of the paper with the exception of the signifi- 
cant wave height which was 10 ft instead of 14 ft 
(the significant wave height appropriate for Fig. 
6). The Rayleigh distribution for this case 
showed that the pitch angle having the highest fre- 
quency of occurrence was about 1.9 deg. If 
linearity of response can be assumed, the most 
probable pitch angle when the ship is running in 
waves whose significant height is 14 ft would be 
2.7 deg. This appears to be in close agreement 
with the results of the author’s Fig. 6. 

One further remark may be of some interest. 
Data obtained on the pitching motions of a trans- 
port ship also were fitted by the Rayleigh dis- 
tribution. The data were collected at %-hr in- 
tervals over a 36-hr period. The significant wave 
height varied from approximately 5 to 15 ft and 
the wave period varied from 5 to 12 sec. The 
wave direction during this period changed by 3.5 
deg. The average speed of the ship was 17 knots. 
While, according to the author a log-normal dis- 
tribution should more accurately describe the 
data, the Rayleigh distribution was found to pro- 
vide a very reasonable fit in this case also. 

Finally, one might remember that while the 
author’s method does seem to offer a convenient 
method for organizing a large amount of data, 
the power-spectrum technique, while more labori- 
ous, appears to offer a solution of even greater 
practical significance. 


Pror. W. J. PIERSON, JR.,”” Visitor: This very 
fine paper will certainly provide the naval archi- 
tect with important design criteria in his work. 
The words of warning in the author’s discussion 
are particularly apt, and they should be heeded! 
The author of these comments agrees with a 
major portion of the results of the paper and be- 
lieves that quite a few years will elapse before 
they are superseded by more refined analysis. 

A few points need to be made, however. 
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The author states, ‘‘the experimental approach 
taken in this paper should complement the theo- 
retical work of St. Denis and Pierson, although the 
present study was developed independently of 
theirs and is not limited by the assumptions re- 
quired by their mathematical analysis.””’ In my 
opinion, all of his results are strictly limited by the 
assumptions of the mathematical analysis em- 
ployed by Mr. St. Denis and the writer. In par- 
ticular our assumption that the sea surface could 
be represented by a stationary Gaussian process 
and that the response of the ship to the waves is 
linear is at the very root foundation of the Ray- 
leigh distributions studied in this paper. It also 
follows that these assumptions are at the founda- 
tion of the log-normal distribution. 


It is interesting to note that the Rayleigh distri- 
bution and the log-normal distribution fit so well. 
However since a stationary Gaussian process is a 
linear process, since waves are definitely nonlinear, 
and since extreme ship motions are nonlinear, it 
follows that neither the log-normal nor the Ray- 
leigh distribution will truly describe conditions at 
extremely high values. That no such discrep- 
ancy was detected in this work is probably due 
to observational error and sample size, but never- 
theless it must exist. 


Asan example, consider a sea generated by a 30- 
knot wind blowing for an infinite length of time 
over an infinite fetch. A-wave record from such a 
sea would have a definite E value. If one waited 
long enough, according to the extreme-value 
theory of Longuet-Higgins, a wave higher than 
any preassigned arbitrarily high value would pass. 
However, the sea surface is nonlinear, and the 
spectrum appears to have no energy below a cer- 
tain frequency. Therefore waves in excess of a 
certain height cannot exist because they would 
break. Our knowledge of the process is not suffi- 
ciently precise to tell us whether the limiting non- 
linear wave is given by the highest wave out of 
1000, 10,000, 100,000 or 1,000,000 according to 
the statistics of Longuet-Higgins. Similar re- 
marks are true about any other quantity with a 
Rayleigh distribution in this paper. 

Similarly there is an upper bound to the wind 
speeds possible over the ocean, to the duration of 
the wind and to the length of the fetch over which 
it blows. And there is a highest possible wave 
due to nonlinear considerations for these condi- 
tions. The log-normal distribution may well be 
too high for higher values for the foregoing 
reasons. 


The Rayleigh distribution fits a great deal of 


these data quite well. Since the “Chi Square”’ 
test gives the probability that a sample truly 
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drawn at random from a population actually with 
the probability density function being tested could 
differ by chance by as much as it does from the 
theoretical distribution, one can argue that charac- 
terizing the data as in Table 3 by terms such as 
“fair fit’ or ‘‘very good fit” is not justifiable. One 
can only state that the hypothesis that the data 
were drawn from the theoretical population being 
tested cannot be rejected at the 5 per cent signifi- 
cance level. It is odd that in one hundred such 
tests of ship motion or hull stresses no discrep- 
ancies at the 5 or 1 per cent level were found. 
One should expect to reject the hypothesis at the 
5 per cent level about 5 times in 100 trials and 
once at the 1 per cent level (despite the fact that 
it is true). 

Recently Cartwright and Rydill?* have found 
that the Rayleigh distribution did not fit certain 
ship-motion records and wave records which they 
were studying due to the fact that the spectrum 
was not anarrow band. They give a more refined 
formula which fits the data better. 


Dr. J. H. CHapwick, Associate Member: It has 
been said often that the main object of scientific 
research is to discover and explain regularities in 
the data. In this sense the paper is certainly a 
contribution to marine science as well as to marine 
engineering—for what could be more regular than 
a straight line? This paper will be, I am sure, 
another important landmark on the trail recently 
blazed by St. Denis and Pierson in the paper they 
presented before this Society in 1953. Where 
St. Denis and Pierson emphasized theory and fre- 
quency distribution (power spectrum), the present 
author emphasizes fact and amplitude distribu- 
tion. Taken singly or together, these papers 
offer us an infinitely more realistic view of ocean 
waves and ship motions than we have had hereto- 
fore. 

The two papers are closely related and comple- 
ment each other in a natural way. I have to 
violently disagree with my friend Dr. Pierson, 
however, when he argues that Jasper’s results are 
based on his linear model assumptions. The 
author’s results are largely factual, and it is a 
fundamental principle of science that theories are 
based on assumptions, while facts are based on 
measurements. In support of this view may I 
quote the epigram that appeared at the head of 
the original St. MDenis-Pierson paper (15): 
“Mathematics can never tell us what is; only 
what would be if.’’ (Poincare). It is true that 
some of the author’s results can be explained 
(i.e., predicted) in terms of the linear model 
assumptions; which is to say that his facts so far 
tend to confirm the. St. Denis-Pierson theory. 
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Nevertheless, theories are not facts, and never 
the twain shall meet. 

Both of these papers have as a central theme, 
the notion that waves and wave-induced ship mo- 
tions are stochastic processes; i.e., essentially 
random in nature. Actually, this is something 
that marine science has recognized implicitly for a 
long, long time, but we have stood shivering on the 
verge of an idea, unwilling to plunge into the seem- 
ingly infinite mathematical complexities of the 
statistical approach, until Messrs. St. Denis and 
Pierson came along to push us in. A surprising 
number of us survived this immersion. To the 
few that have felt themselves going down for the 
third time, the author now throws a life-ring. 

Mr. St. Denis, in another paper on ship mo- 
tions, once quoted Lord Rayleigh as saying: ‘“The 
basic law of the seaway is the apparent lack of any 
law.” It is pleasantly ironic that one basic law 
of the seaway, discovered during the evolution 
of these new theories, is that peak amplitudes over 
the short term tend to fall on a Rayleigh distribu- 
tion, a fact not known in Rayleigh’s lifetime. This 
little joke of fate tells the whole story. The basic 
law of the seaway is randomness, but randomness 
itself obeys statistical laws, some of them as 
simple as straight lines. 

The author stresses engineering applications, 
and the point is well taken. Since the present 
models are vastly more realistic than the old single- 
sine-wave models they also have a much greater 
potential value to the engineer. The amplitude- 
distribution data presented in the paper are par- 
ticularly important with respect to the economic 
side of design. 

In the collection of his data, the author has used 
a technique the importance of which cannot be 
overemphasized; namely, the collection of data at 
predetermined times. This technique provides 
objective data with the minimum intervention of 
theory, and hence provides the best possible test 
of theory. In its application, this technique has 
certain consequences: It tends to call for a long- 
term, low-level effort; it requires automatic re- 
corders; and it uses a minimum of shipboard per- 
sonnel. At the other extreme we have the grand 
expedition approach, using short-term, high-level 
efforts, subjective recording and a maximum of 
shipboard personnel. History has shown that 
the expedition approach is an extraordinarily in- 
efficient method of collecting data about wave- 
induced ship motions. The San Francisco and 
Ocean Vulcan expeditions are cases in point, par- 
ticularly the former. It is to be hoped that in the 
near future we will see the same type of recording 
technique used by the author for amplitudes ap- 
plied also to power spectra and cross spectra. 


Mr.R.P. Erickson,” Visitor: The excellent 
statistical description of ship data, discussed in 
the paper, has as a major application the establish- 
ment of design criteria. Because the data are 
empirical they fit ideally only the ship on which 
the data were gathered and must be extrapolated 
for new ship designs. If these new designs are 
not distinctly different from the tested design it 
follows that the extrapolation would be satisfac- 
torily accurate. The purpose of this discussion is 
to mention a problem where previous similar ex- 
perience is not available and thereby point out a 
need for the statistical description of derived, as 
well as empirical, data. 

This problem arises in the specification of en- 

_ vironmental conditions encountered during the 
shipboard handling and launching of large ballis- 
tic missiles. The structural loads imposed on the 
missile by the motion of the ship is one of a number 
of similar problems that represents a critical de- 
sign criteria. Minimum structural weight is abso- 
lutely essential to the effective utilization of the 
missile as a weapon. Minimum weight implies 
minimum acceptable structural strength. The 
emphasis on weight reduction is so strong that 
missile designers are prepared to accept a statis- 
tical risk of structural failure if a known statistical 
risk can be defined. 

The establishment of a statistical design criteria 
for ship-launched ballistic missiles must be based 
upon an intelligent specification of the maximum 
operational sea conditions in which the weapon 
must be utilized. Then attempts must be made 
to determine the statistical description of the sea, 
or preferably of the ship motions, encountered 
under this assumed steady environmental condi- 
tion. Thus far the attempt to define a design ap- 
proach can benefit significantly from the work of 
the author. Now, however, because of a lack of 
previous experience, additional information as to 
missile structural loads must be derived by en- 
gineering analysis and the statistical properties of 
the derived information must be determined. 

To date, with the co-operation of the Navy, 
especially Bu Ships, the Hydrographic Office, and 
the David Taylor Model Basin, we have been able 
to define statistical wave-height and wave-length 
distribution based on constant underlying param- 
eters. We have been forced to extend this in- 
formation to the statistical definition of missile 
loads by relatively crude graphical methods. 
This technique can be improved by additional ex- 
perimental data on ship response which is now 
being gathered. The eventual goal of this pro- 
gram is to be able to describe derived data in a 
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convenient statistical form similar to that sub- 
stantiated by the author for experimental data. 


Mr. Harry Press,”?? Visitor: The author is to 
be congratulated in making a valuable contribu- 
tion to the problems of ship behavior in disturbed 
seas by his collection and analysis of statistical 
data on ocean waves and ship responses. The 
paper is of particular interest to me in that it 
closely parallels some of our recent studies on the 
analogous problem of the airplane behavior in 
rough air published in the National Advisory Com- 
mittee for Aeronautics Technical Notes 3362 and 
3540. 

The paper considers two major points: The first 
point is concerned with the demonstration of the 
applicability of the Rayleigh distribution for de- 
scribing the probability distribution of ocean 
wave heights and ship responses for steady ocean 
conditions. The second point is concerned with 
the probability distribution of wave heights and 
such ship responses as disturbed motions and 
structural stresses under broader operating con- 
ditions, presumably over-all ship operations. The 
author claims that such distributions are well 
represented by a log-normal probability distribu- 
tion. In both cases, the evidence presented is al- 
most entirely empirical and the conclusions 
reached based on statistical tests on the data of 
these hypotheses. 

A few comments on the two major points of the 
paper appear warranted. In regard to the Ray- 
leigh distributions for steady ocean conditions, 
the statistical evidence presented in most cases 
seems to be quite convincing. It is unfortunate 
that the author didn’t more closely tie the present 
statistical evidence to the closely related recent 
work on the response of linear systems to Gaussian 
random processes. Results obtained in the theory 
of stationary processes indicate that the distri- 
bution of peaks for a Gaussian process having a 
narrow power spectrum should tend to adhere to a 
Rayleigh distribution (Rice). The work of Pier- 
son and St. Denis has served to demonstrate that, 
to a very good first approximation, the disturbed 
sea under given conditions may be considered a 
stationary Gaussian random process. Conse- 
quently, if the ocean wave power spectrum is 
narrow and peaked, which appears to be the case 
frequently, the distribution of ocean wave heights, 
either peak to trough or mean value to peak, may 
well be expected to approximate a Rayleigh dis- 
tribution. This approximation would tend to 
deteriorate as the spectrum broadens. 

In regard to the distribution of peaks of the ship 
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response such as pitch, structural stresses, etc., 
the distribution of peaks may be expected to de- 
pend upon the characteristics of the ship’s fre- 
quency-response functions assuming that the ship 
isa simple linear system. As long as the ship re- 
sponse involves a single degree of freedom that is 
relatively lightly damped (and thus yielding a 
narrow and peaked frequency-response function), 
the Rayleigh distribution may still be expected 
to apply. However, if the ship response is more 
complicated involving several significant degrees 
of freedom, the resulting power spectrum of the re- 
sponse may no longer contain a simple peak but 
may contain multiple peaks. For such cases, the 
Rayleigh distribution would not appear to be ade- 
quate, and the more general distribution given by 
Rice would have to be used. 

The writer is not familiar with how serious the 
foregoing objections actually may be in practical 
applications. The experimental data presented 
do, however, show sufficient departure from the 
assumed hypotheses to warrant a more detailed 
consideration of the foregoing types of limitations. 

In regard to the second major point, the ap- 
plicability of the log-normal distribution, we are 
here concerned with a distribution of peaks for the 
summation of a large number of stationary Gaus- 
sian processes. There does not appear to be any 
theoretical a priori basis for the choice of a distri- 
bution function, and one is forced here to rely more 
heavily on the purely empirical observations. 

The empirical data show generally good ad- 
herence to the assumed log-normal distribution 
function, but again the data appear to show de- 
partures from log-normal distributions that sug- 
gest the possibility of the need for a more detailed 
consideration. It appears clear that the appro- 
priate distribution may well be expected to vary 
depending upon the detailed nature of the opera- 
tions being considered and the climatology of the 
ocean waves. Inananalogous problem in airplane 
studies, no single distribution was found to apply 
in general to gust velocities encountered by air- 
planes and to airplane response distributions under 
full operating conditions. 

One further point is worth noting. The author 
has chosen to describe the ocean characteristics in 
terms of the probability distribution of the ocean 
wave heights. Although this is a useful descrip- 
tion for many purposes, it is not adequate in re- 
gard to the problem of determining the ship re- 
sponse since the response of a ship will depend 
upon not only the individual wave heights but 
also upon the preceding waves or more concisely 
upon the power spectrum of the ocean wave 
height. Thus, a more fundamental and useful de- 
scription of disturbed seas would appear neces- 


sary in the long run if adequate predictions of ship 
behavior in operations are to be achieved. 

In concluding, I should like to reiterate the 
valuable nature of the present paper. It cer- 
tainly provides some reasonable and useful defini- 
tions of the ocean wave field and the character of 
the ship response. The principal criticisms in the 
foregoing are concerned with the absence of ex- 
pressed limitations to the hypotheses claimed and 
the suggestion is that refinements to the simple 
picture presented may often become necessary. 


Compr. R. L. Brooxs, USN, Associate Mem- 
ber: This paper provides a working tool of great 
value to the naval architect for obtaining, in- 
terpreting, and using full-scale data. By means of 
the distribution pattern and with a minimum of 
observations, many seaworthiness problems can 
be clarified in terms of actual results obtained at 
sea on a properly chosen hull. 

That the need for reliable data respecting ship 
responses to the sea is a matter of great impor- 
tance has become better appreciated, as the solu- 
tion of each new seaworthiness problem was found 
to depend on full-scale data which were seldom 
readily available in the form or ranges which were 
needed. Examples of critical problem areas in 
which progress no longer needs to be hampered by 
the lack of full-scale data, now that this distribu- 
tion pattern technique is available, include: (a) 
The determination of how much ships must slow in 
heavy seas, (b) the checking of assumptions upon 
which are based theories for predicting ship 
motions, (c) the improving of present strength 
design methods, and (d) the development of ship- 
board computers for selecting optimum courses 
and speeds for a variety of tactical situations, 
as for example fueling of destroyers in moderate 
to heavy seas. 

While the need for reliable data has long been 
acknowledged, the difficulties which stood in the 
way of obtaining it at sea, have been so serious 
as to hold up solutions to many of the important 
problems for several years. 

Consideration of a few of these difficulties helps 
to show that the present paper represents a task 
of large magnitude, because each difficulty was a 
major problem in its own right. The following 
examples are mentioned: (a) The development 
of specialized recording instruments which would 
require a minimum of attention, (b) the planning 
of how the many variables which define the 
environment and the ships’ responses should be 
handled, (c) the detailed co-ordination of in- 
strument installation and operation, so as to 
harmonize with ship schedules and internal 
routine, (d) the difficulty of defining each sea 


52 


state in a manner which would have consistent 
meaning to all who would be involved in the work 
or would use the results, and (e) the presentation 
of data in a manner which would lend itself to a 
maximum of usefulness. 

The work involved in this paper required sev- 
eral years of patient, painstaking, and persistent 
effort; and called for great vision and faith on the 
part of the author. He not only has succeeded in 
developing the necessary instruments and meth- 
ods, but also has already used them on many 
ships, and has applied the tools of statistics and 
oceanography in collecting and presenting the 
results. These much-needed accomplishments 
are deserving of great credit to the author. Much 
credit also is due to the several organizations 
which nurtured this work, in particular to the 
Society, the Navy, the Model Basin, and the 
Coast Guard. 

With the advent of nuclear power for driving 
surface ships, it is well that this paper gives us the 
tools for defining and correcting the weaknesses 
of hull designs, before they become handicaps 
under nuclear speed requirements. 

Admiral Rickover has stressed this point in a 
recent issue of Naval Institute Proceedings in the 
following manner: ‘With nuclear power, the 
Navy can go where it wants, when it wants, at 


top speed without worrying about) fuel economy. 


However, this highly desirable goal cannot be 
attained unless the hull is also adequate as to its 
seaworthiness characteristics. Thus, nuclear 
power undoubtedly will bring greater emphasis 
on these three questions: (a2) How much of the 
speed potential must be given up in medium and 
heavy seas? (0) If the hull is inadequate, in what 
way must its design be changed? (c) To what 
extent do physiological and psychological factors 
influence the decision to slow ships in heavy 
seas? The answers to these questions are im- 
portant for conventionally powered ships and be- 
come even more so for those with nuclear drive. 

It is hoped that work in this field will be carried 
forward vigorously; and particularly that the 
problems of determining distribution patterns by 
model testing, and of determining fully the in- 
fluence of alternating stresses on hull structures 
and the limitations as to heavy sea operation, 
whether from stresses or motions, will be com- 
pleted within the next few years. 


Dr. JuLius Liesiein,® Visitor: It is under- 
stood that some comment has been raised con- 
cerning a statement in the paper to the effect that 
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statistical tests of significance showed no signifi- 
cant deviations in any of the more than one 
hundred distributions tested. The argument 
raised was that even with a significance level of 
say, 95 per cent, about 5 per cent of the tests 
should show significance even when none is pres- 
ent. 

First, it may be noted that the statement of no- 
significance cases may have been based upon 
visual inspection in conjunction with calculation 
of some values of “Chi-Square.’’ This could 
allow a few cases to be actually significant with- 
out appearing so to the eye. 

Secondly, for the practical engineering pur- 
poses intended, the postulated Rayleigh distribu- 
tion seems amply supported by the data and past 
experience. In fact, the quality of statistical 
treatment is above that usually to be found in in- 
vestigations of this character, and the issue of 
combining numerous statistical tests is rarely, if 
ever, raised—perhaps because the data seldom 
allow calculations to be carried that far. More- 
over, the too few significance cases might not be a 
rare occurrence. It may be pointed out that ex- 
tensive undertakings are often carried out on the 
basis of conclusions obtained without detailed in- 
vestigation of the results of significance tests. 


Pror. E. J. GuMBEL,*! Visitor: Being a statis- 
tician, I want to talk only about the statistical 
aspects of this paper. The graphical way he 
used the Rayleigh distribution is very able, be- 
cause he reduces it to the exponential distribution. 
This simplifies considerably the plotting, because 
a semi-logarithmic paper can be used. 

The main point of his contribution is the state- 
ment that “distribution of wave heights,’ and 
many connected phenomena, ‘“‘may well be ap- 
proximated by the Rayleigh distributions, pro- 
vided that environmental conditions are steady 
and that the same variables will be reproduced by 
a logarithmically normal distribution if the under- 
lying environmental conditions are allowed to 
vary over a wide spread.” 

The first statement has a theoretical basis; 
while the second statement is purely empirical. 
This leads to the question how to prove that a con- 
volution of Rayleigh distributions leads to a 
logarithmic normal distribution. 

The logarithmic normal distribution is not 
(as the author and bad textbooks claim) the nor- 
mal distribution of the logarithm. In reality, 
the probability function is an integral and the dis- 
tributions are transformed as functions under an 
integral. The incorrect statement, however, has 
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no influence on the merit of the paper or on the 
results, because only the probability function is 
used. 

If we define the significant wave heights as being 
the mean third highest wave, then the distribu- 
tion of the third highest value is a modification of 
the distribution of the largest value. The proba- 
bility points for the distribution of the third 
highest values are given in the Bureau of Stand- 
ards Tables on extreme values. 

The distribution of the largest value and of the 
third largest value have some geometric similarity 
with the logarithmic normal distribution. But 
in our case, it is more natural to use the distribution 
of the third largest value for the significant wave 
heights. The means of the third highest value 
should be normally distributed, just as the means 
of the largest value are normally distributed. 

I agree with the final statement: If the stresses 
are logarithmically normally distributed, then the 
largest stress has the first asymptotic distribution 
of the largest value. These remarks are given 
only from the statistical standpoint and should 
not be considered as a critique of the intrinsic 
technical value of this paper. 


Dr. JASPER: The wide interest shown in this 
paper is greatly appreciated and the author wishes 
to thank the discussers for their contributions 
which will add considerably to the value of the 
paper. Several discussers (Professors Korvin- 
Kroukovsky and Pierson, Dr. Szebehely, Messrs. 
Press and Marks) indicate the desirability of pro- 
viding a theoretical foundation for the proposed 
distribution patterns. The author fully agrees 
that this would be desirable but is neither in a 
position to do this adequately nor to evaluate 
available theories of ship response to waves, such 
as have been proposed by St. Denis and Pierson, 
on the basis of the measurements reported here. 
Therefore, there would appear to be little reason 
to elaborate on these well-known theories at this 
time. 

The distributions given here are based on 
physical measurements, i.e., facts; the statistical 
theory is utilized to describe these facts concisely 
and to test if there are significant differences be- 
tween the assumed statistical distribution and the 
factual data. Therefore the proposed distribu- 
tions are derived from experiment alone, that is 
they are empirical and therefore are not limited 
by prior assumptions such as always are the basis 
of theories. This is contrary to the opinions ex- 
pressed by some discussers. Dr. Chadwick's 
comments on this matter are very much to the 
point and are fully concurred in. 


It is not precluded, as pointed out by several 
discussers, that mathematical distributions may 
be found which fit the data better than those sug- 
gested here, though no doubt at the expense 
of additional complexity. Undoubtedly the true 
physical distribution is much more complex than 
any discussed here but this does not prejudice 
the practical utility of a much simpler distribu- 
tion. 

It has been stated by Mr. Lewis and others 
that the Rayleigh distribution applies only to 
time-stationary processes with narrow spectra 
and that therefore the author’s work indicates 
that most sea and ship-response spectra are 
“narrow.’’ Objections can be raised to the term 
“narrow” spectrum even if the assumption of a 
stationary time series process is accepted, al- 
though the applicability of this or other dis- 
tributions: is obviously not restricted to time- 
stationary processes. A paper by Cartwright and 
Longuet-Higgins*? has demonstrated that, for 
the representative cases studied by them, no sig- 
nificant deviation between the Rayleigh distribu- 
tion and test data was apparent for spectra with a 
“root mean square width” e of less than 0.48, 
where « = 0 corresponds to the exact Rayleigh 
distribution and e = 1 to the normal distributions. 
Although they did find significant deviations from 
the Rayleigh distribution for some ocean-wave 
records no such deviations were indicated for 
the ship’s response to the waves. We might then 
define ‘‘narrow” spectra as those for which e < 
0.50. 

Several discussers (Marks and Pierson) have 
mentioned the recent paper by Cartwright and 
Rydill?? which makes some comparisons be- 
tween distribution patterns obtained by theory 
and by sea tests. The statistical function used by 
them to represent the short-term distribution of 
ship motions is more general than the Rayleigh 
pattern. They applied it to the distribution of 
the maximum deviations from a mean value 
rather than to the peak-to-peak variations dealt 
with in the present paper. It should be noted 
that these two variables in general would not be 
accommodated by thesamefunction. The authors 
claimed that there were no significant differences 
between the theoretical distribution and the ex- 
perimental histograms but did not give the basis 
for this claim. Personal correspondence with 
the authors indicates that the question of good- 
ness of fit had been somewhat neglected and that 
there were indeed significant departures from the 
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assumed distribution. However they felt that no 
valuable (statistical) conclusions could have been 
drawn from their results since they were not really 
in a form suitable for vigorous significance tests. 
Nevertheless, the short-term distribution used by 
them?*2 probably would be theoretically more 
accurate than the Rayleigh distribution. 

Professors Pierson and Korvin-Kroukovsky 
have indicated that the fit of the proposed dis- 
tribution to the data is ““‘unusually’’ good whereas 
Messrs. Lewis and Press note occasional appreci- 
able departures of the data from the theoretical 
distributions. The fit is really not “too good to be 
true.” Dr. Lieblein’s discussion bears directly on 
this question. 

It should have been stated in the paper that 
not every one of the hundred or so distributions 
was actually subjected to statistical tests. Rather, 
all these distributions were examined visually and 
then several typical distributions as well as sev- 
eral relatively poorly fitting distributions were 
selected for statistical testing. This seemed to be 
a reasonable procedure. Also the confidence 
limits fitted to the cumulative distribution are 
applicable individually to each experimental frac- 
tile but they do not by their nature extend to very 
small or very large fractiles, unless very large sam- 
ples are used. Thus nothing can be said about 
the fit at these extremes and it is therefore quite 
possible that occasional significant deviations 
could indeed be present at these extremes, with- 
out being detected. Of course from the stand- 
point of practical application all that is desired 
is a fit good enough to permit making reasonable 
estimates. The latter criterion has been satis- 
fied, it is believed, to a good degree for the ship’s 
response to the sea. The author cannot conclude 
that the deviations from the proposed distribu- 
tions, which were noted by Messrs. Lewis and 
Press, are statistically significant. 

Messrs. Lewis, Comstock, and Vasta have em- 
phasized the importance of being able to estimate 
the extreme value of hull girder stress, raising 
the question as to how this could be done ade- 
quately. The long-term distribution would not 
be suitable for this purpose, as stated in the paper. 
The following method is suggested as reasonable: 

Determine £,, for the Rayleigh distribution of 
stress corresponding to the worst combination of 
sea state, ship speed and heading that a ship is 
likely to experience, utilizing either model tests, 
full-scale data, or theory. Estimate the longest 
duration that the ship will operate continuously 
under this combination of steady conditions and 
note the number JV, of variations experienced in 
this length of time. Next, estimate the number 
of times n, that these extreme operating condi- 
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tions (storms) will be experienced during the serv- 
ice life of the ship. Thus (JV) stress variations 
are expected corresponding to a basic distribu- 
tion defined by E,. Then, utilizing Formula 
[2] of the paper, we can estimate the fraction f 
of all similar ships which, on the average, will ex- 
ceed om over their service life. Thus f represents 
the risk of exceeding om. Since for the Rayleigh 
distribution 


[1 — P(on)] = e7*n'/En 
we have from Formula [2] 
1 — f = exp[—e~] 


On" 1 N 
y= —= => Oe (n ) 


Table 1 of Bureau of Standards bulletin ‘“Ap- 
plied Math Series No. 22” tabulates the value of 
exp [—e-’] as a function of y. Thus for a 
specified risk f of exceeding the stress om one may 
look up the corresponding value of y and then 
solve for the desired stress from the relation 


Om? = Emly + loge(Nn) ] 


Another way of making the estimate is to 
utilize the approximation that if f; is the proba- 
bility of exceeding om in a single storm then 
nf; is the probability of exceeding om in m storms. 
That is f = fi. The formulas given here thus 
far are valid only if f is much smaller than unity 
which will be true if we are looking for extreme 
values. The exact expression for f is 


foe loa@oipp 


As an example consider an aircraft carrier with 
JB 16, N = 1000, x = 100, and take f = 
0.001; i.e., one chance in a thousand. Then y 
7.0 and 


om = [16(7.0 + 11.5)]!/? = 17.2 kips/sq in. 


The method of Hazen and Nims referred to by 
Mr. Lewis is likely to result in an overestimate of 
the stress but otherwise appears to offer an at- 
tractive solution. These approaches may be 
utilized to establish an upper limit to the long- 
term distribution also. 

The short-term Rayleigh distributions fitted by 
Mr. Lewis to experimental model test data of wave 
height, ship motions and bending moments (Figs. 
24, 25) do show encouraging agreement. It is the 
author’s expectation that, once agreement be- 
tween model tests and full-scale tests has been 
established, model tests will be used to define the 
distribution patterns. Mr. Lewis’ fear that an 


underestimate of the extreme wave-induced hull- 
girder stress would “inevitably doom many ships 
to the bottom”’ is believed to be rather on the 
pessimistic side. There are many other stresses 
present besides the wave-induced stress and not 
all will be likely to experience their maximum 
value at the same time. 

The comments of Mr. Darbyshire are par- 
ticularly welcome as he has pioneered in the 
measurement of ocean waves with a shipborne 
wave recorder. His data provide the first ex- 
tensive series of measurements of ocean-wave 
heights. It is very gratifying to note that Mr. 
Darbyshire finds”! that the proposed log-normal 
long-term distribution fits the ocean data satis- 
factorily. It might be well to note at this point 
that the author did not study the short-term dis- 
tribution of wave heights and made no attempt 
to show that their distribution was of the Ray- 
leigh type; rather he stated that “it has been 
accepted generally that the heights of ocean waves 
passing any given point follow the Rayleigh dis- 
tribution. ..’’ Recent studies indicate that of all 
the variables discussed here, the short-term pat- 
tern of the height of the wave surface is least well 
represented by the Rayleigh distribution. 

Both Mr. Comstock and Mr. Vasta express the 
opinion that endurance strength is not a problem 
in ship structures, inferring that only the maxi- 
mum total stress need be considered. It is not 
believed that there is a sound basis for this con- 
clusion. A study by Dr. Weck** has shown that 
fatigue failures of typical welded plate-stiffener 
combinations may occur at very low nominal 
stresses (+ 2 tons/sq in.) at the point of failure. 
Such cracks, when noted, probably would be re- 
paired without coming to the attention of de- 
signers; they could on occasion form the nucleus 
of a major fracture in the presence of large tem- 
perature-induced, dead-load, or seaway stresses. 

The discussions of Weck’s paper force one to 
one of two conclusions, quoting Dr. Weck: 
“Hither that fatigue failures do not occur in ships, 
which seems to be contradicted by the majority 
of those who have contributed to the discussion, 
or that they are not recognized as such...” It 
also may be timely to quote the following from his 
paper: “It is difficult to see why it should be 
profitable to distinguish between the main hull 
structure of a ship and the rest when considering 
the incidence of fatigue. This distinction may 
well be convenient for purposes of design, but it 
does not correspond to reality.” 

“A study of ship failures sometimes referred to 
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as providing a basis for concluding that fatigue 
failures are not important is the “Final Report 
of a Board of Investigation” dated 15 July 1946, 
which was made for the Secretary of the Navy to 
study the high incidence of failures in welded ships. 
The service data given in that report are heavily 
weighted in favor of failures incident to the brittle- 
fracture problem because the ships were es- 
pecially subject to this ‘‘disease’’ and thus could be 
expected to “die” thereof before an endurance 
failure could develop. What is needed is a con- 
sistent and thorough system of reporting struc- 
tural failures and ‘‘defects’’ even though they may 
not appear dangerous. The author did not in- 
tend to convey the impression that “‘fatigue’’ is a 
major structural problem in ship design—he does 
not know—however, this problem should be 
studied together with the other possibilities of 
causing failure, without prejudice. 

The author is not in a position to answer Mr. 
Vasta’s question as to whether hull-girder struc- 
tures are designed too conservatively. as he has no 
sound basis for a definite conclusion. However 
the author “‘feels’’ that the sectional moment of 
inertia is probably more than adequate in many 
ships. On the other hand the constructional de- 
tails with which the structural material is put 
together may not always be adequate and further- 
more the local strength, such as that of bottom 
structures subjected to slamming, is on occasion 
insufficient. It is really not realistic to separate 
local strength from hull-girder strength. The 
author fully appreciates Mr. Vasta’s other com- 
ments and is in general agreement with them. 

Mr. Marks is increasing the scope of the paper 
when he states that it attempts to solve the ship- 
performance problem. This is not so. The ob- 
ject, as has been stated repeatedly, is to establish 
that the long and short-term distribution of ship 
motion and stresses approximate to the Rayleigh 
and log-normal distributions, respectively, and to 
give illustrative engineering applications of this 
result. True, it was hoped that a knowledge of 
the form of the distribution would simplify the 
solution of the ship-performance problem. As 
stated in the summary, it is expected that this so- 
lution will take the form of model tests in waves 
which will be utilized to make predictions for par- 
ticular ships in a particular service analogous to 
the procedure now used to determine power re- 
quirements in smooth water by model tests. 

Professor Korvin-Kroukovsky feels that in the 
paper there is a lack of reference to the physical 
nature of the situation. The physical nature is of 
course reflected in the parameter E (the mean 
square value of the physical motion) which re- 
lates the mathematical distribution to the “physi- 


cal situation.”” A paper by Williams was re- 
ferred to as an illustration of a physical basis for 
the Rayleigh distribution. The mathematical 
expression given by Williams for a ‘fluctuating 
oscillation” together with the stipulated assump- 
tions made by him does give rise to the Rayleigh 
distribution. But what is the advantage of 
carrying out such a derivation? If the mathe- 
matical model did describe the physical situation 
correctly and generally, this would indeed prove 
the Rayleigh distribution to be the applicable dis- 
tribution function. However, there appears to 
be no reason to believe that this mathematical 
model is the correct one, and if it is not, then the 
mathematical development is superfluous. 

It was the intention in this study to work di- 
rectly with service motions and stresses of ships 
thus arriving at results which do not depend upon 
the validity of a priori assumptions. Thus there 
has been no attempt made here to relate cause 
and effect. If a theory relating the sea to the 
ship’s response should be proven practical to 
apply, then distribution functions would still have 
to be established. Thus it seemed prudent to 
start by making an attempt to establish these 
probability distributions directly since they would 
be of immediate utility, irrespective of the avail- 
ability of any suitable theory on waves and the 
ship response to waves. 

The comments of Professor Pierson have, for the 
most part, already been discussed. The author 
agrees that the proposed distributions cannot be 
valid for very extreme values; this limitation was 
stated in the paper. ‘fhe limitation is however 
not serious in many practical applications. For 
example, steady-state conditions of ship speed, 
heading, and so on, will not generally be main- 
tained long enough to result in predictions which 
are too extreme. 

The remarks of Dr. Chadwick are much ap- 
preciated especially so because he, in a way, 
represents the consumer. He is in a position to 
utilize the information presented here and it is 
pleasing to hear his fresh and forthright com- 
ments on the paper and on some of the points 
raised by others. In fact, Dr. Chadwick has 
effectively answered some of the points raised by 
other discussers. 

Mr. Erickson’s problems are typical of those 
met more and more often in recent years. The 
question of specification of environmental condi- 
tions as to rigid-body motion, vibration, and 
shock may well benefit by consideration of statis- 
tical methods for the collection, presentation, and 
analysis of the data. In some of these areas we 
are not as yet much concerned with a high degree 
‘of accuracy in specifying the environmental con- 
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ditions but rather with assigning a reasonable 
order of magnitude and the probabilities asso- 
ciated with them. 

Dr. Szebehely is amazed “‘by the fact that such 
a large amount of data could be organized accord- 
ing to relatively simple basic laws.’’ He probably 
means to imply that the basic laws are really 
much more complex. The author does not doubt 
that this is true and has not claimed that the 
Rayleigh and log-normal laws are applicable basic 
laws. He merely takes the position that these 
simple laws give results which are good enough to 
be useful and possibly just as good as would be 
obtained if the correct law were known, consider- 
ing the statistical variations together with the 
various unknowns and inaccuracies that enter the 
problem. 

The author is reasonably satisfied to accept the 
good fortune of having simple laws available which 
describe the situation fairly well without knowing 
why this is so. Life and Nature are full of mys- 
teries—if we should refuse to accept their mys- 
teries unless given a scientific explanation, prog- 
ress would be at a snail’s pace. 

I shall attempt to answer several of Dr. Szebe- 
hely’s questions. A number of statistical laws 
were tried before testing the log-normal law. The 
latter was tried because it has been found appli- 
cable in a number of other physical phenomena. 
As far as I can determine there is no definite re- 
lation between the short-term Rayleigh distribu- 
tion and the long-term log-normal distribution on 
a purely mathematical basis unless one assumes 
an a priori distribution of Rayleigh distributions. 
Nature of course provides a distribution of short- 
term distributions and it may eventually be possi- 
ble to establish the connection then. 

The question raised as to the manner of apply- 
ing the results of this paper to the evaluation of 
ship performance could be discussed at length. 
However, briefly the following may be suggested: 

(a) Assuming the model test procedures are 
shown to be valid then model tests may be utilized 
‘to determine E for short-term distribution. The 
parameter F is an index of performance. Thus a 
single number may be a good criterion of per- 
formance rather than a graph or a table of values. 

(6) Long-term distributions of ship response 
may be prepared by synthesis, utilizing typical 
short-term distributions of ship response each 
specified by the parameter E, together with the 
long-term distribution of sea conditions given in 
the paper. This has been done in TMB reports. 

(c) Expected extreme values can be given on 
the basis of the parameter / and the expected 
number of variations. 

(d) Measurements on both model and full 


scale could be made directly in terms of the 
parameter E rather than measuring the variable 
as a function of time. 

Dr. Szebehely’s question regarding the long- 
term distribution applicable to a destroyer shown 
in the paper can be answered as follows: (a) The 
line would be different for different operating con- 
ditions or for a substantially different destroyer. 
It would not be necessary to measure stresses for 
long periods of time to evaluate the relative merits 
of various ship designs. One would need to find 
the parameter E for representative operating con- 
ditions utilizing either model or full-scale tests or 
theory. The long-term distribution could then 
be synthesized as noted before. Running several 
destroyers side by side in the same sea for 30 
min would give a comparison for the test con- 
‘ditions but not for other operating conditions un- 
less recourse can be had to a valid theory. Further- 
more, it would appear to be a rather impractical 
solution. The best approach, at this time, would 
appear to be in the use of model tests to determine 
and compare the desired distributions. 

Dr. Szebehely’s suggestion of determining the 
‘transfer functions by model tests in regular waves 
and then utilizing the energy spectrum of the sea 
to find the behavior of the ship is, the author be- 
lieves, probably the best immediate solution to 
the entire problem. However it requires that the 
method be shown to be applicable over a suffi- 
ciently wide range of conditions. 

It is heartening to see the good checks Mrs. 
Bledsoe was able to provide between the direct 
application of the Rayleigh distribution and the 
autocorrelation technique. If a rapid computer 
such as the UNIVAC is available then the use of 
the autocorrelation techniques may be just as 
quick and certainly more informative than the 
direct method of computing the mean square 
value of the measured variation. However the 
availability of such computing facilities is a 
severe restriction on the general utility of the 
autocorrelation technique. 

Mrs. Bledsoe indicates that significant and 
average heights, and so on, can be obtained from a 
knowledge of the total energy of a power spectrum. 
This is not believed to be generally true; it is of 
course true for a “‘narrow’’ spectrum. 

The physical example of a distribution problem 
of pitch motion was cited in which the significant 
wave heights varied over.a range of from 5-15 ft 
and the Rayleigh distributions still gave an ade- 
quate representation. This is of course interest- 
ing and significant but it does not in any way con- 
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flict with the applicability of the log-normal dis- 
tribution for representation of long-term dis- 
tribution. The variation of the environmental 
conditions simply did not cover a wide enough 
range. This illustrates the fact that although a 
narrow spectrum gives rise to a Rayleigh dis- 
tribution the reverse is not necessarily true. 
This point seems not to have been given emphasis 
thus far and is pertinent to the comments of 
Professor Pierson. 

The comments of Commander Brooks as well 
as the considerable personal contributions made 
by him to this program during his tour of duty at 
TMB are deeply appreciated by the author. 
Commander Brooks has rightly called attention to 
the possibility that with the advent of nuclear 
power it may occur that full utilization of the 
available speed potential will be limited by the 
lack of adequate seaworthiness characteristics. 
This may well become a major problem facing the 
naval architect in the not too distant future. 

Dr. Gumbel is a well-known authority on ex- 
treme value theory and the author is pleased to 
receive the approval of so well-qualified a statis- 
tician on his approach. The author regrets he 
did not appreciate the fine distinction between the 
logarithmic normal distribution and the normal 
distribution of the logarithm. For the record, 
the definition from Cramer (reference 6) will be 
given: If log x is normal (u,c) then the variable x 
itself has the frequency function 
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This distribution is called the logarithmic normal 
distribution. Thus the same physical situation is 
equally well described by the normal distribution 
of the logarithm of x and by the logarithmic nor- 
mal distribution of x. The only difference is in 
the mathematical expression utilized; i.e., seman- 
tics. In order to keep to the usage of statis- 
ticians the author Has in the paper now given the 
expression for the probability density of x rather 
than the density of log x. 

Dr. Lieblein’s comments are greatly appreci- 
ated because, coming from a statistician, they 
should help engineers evaluate the adequacy of 
the statistical treatment. 

In closing, the author wishes to express his sin- 
cere thanks to all for their interesting discussions. 
It is hoped that some of the questions raised have 
been answered. Any further questions that may 
arise will receive the serious attention of the 
author. 
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